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ABSTRACT 

This  study  was  conducted  to  establish  the  performance  parameters 
of  a  radial  inflow,  dual  discharge  turbine  and  to  determine  the  effect 
of  variations  in  axial  clearance  on  these  parameters.  The  representa- 
tive stream  surface  is  taken  at  the  outer  discharge  radius  instead  of 
at  a  computed  mass-average  discharge  radius,  as  was  done  previously. 
This  technique  results  in  considerably  simplified  computations  and  in 
better  correlation  of  the  rotor  loss  parameters. 

Tests  were  conducted  at  axial  clearances  from  0.015  to  0.081  inches 
and  at  total- to-static  pressure  ratios  from  1.2  to  1.7  for  each  clear- 
ance.  The  test  installation  is  located  at  the  Turbo-Propulsion  Labora- 
tory of  the  Naval  Postgraduate  School,  Monterey,  California. 
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SECTION  I 


INTRODUCTION 


As  requirements  for  smaller  more  powerful  turbines  increase,  to  keep 
pace  with  the  needs  of  the  rapidly  expanding  aerospace  industry,  more  and 
more  attention  is  given  to  radial  turbines.  Since  the  radial  turbine  can 
extract  more  energy  per  stage  than  an  axial  turbine,  it  is  more  useful  in 
applications  where  space  limitations  preclude  the  use  of  a  staged  axial 
turbine.  Currently,  it  is  therefore  used  extensively  for  aircraft  elec- 
tric power  generation,  in  aircraft  cockpit  and  equipment  cooling  systems, 
in  cryogenic  systems,  and  in  missile  auxiliary  power  systems. 

As  a  consequence  of  the  increasingly  widespread  use  of  the  radial 
turbine,  the  establishment  of  performance  parameters  as  an  aid  to  the 
design  of  these  turbines  is  of  fundamental  importance.   It  is  the  objec- 
tive of  this  report  to  establish  some  of  the  performance  parameters  for 
a  particular  radial  inflow,  dual  discharge  turbine  using  air  as  an 
operating  fluid  and  relating  these  parameters  in  such  a  way  that  they  may 
be  applied  to  geometrically  similar  turbines  operating  with  different 
fluids  and  in  different  environments. 

Although  similar  studies  have  previously  been  conducted  on  this 
machine,  this  investigation  differs  in  two  important  aspects.   First, 
several  significant  and  necessary  changes  in  the  internal  dimensions  of 
the  turbine  were  made,  thus  resulting  in  a  somewhat  different  machine 
than  that  previously  tested.  Also,  the  method  of  analysis  used  is  a 
considerably  simplified  one,  which  yields  more  significant  results  with 
less  computational  effort  than  those  methods  previously  employed  by  re- 
ferring all  losses  to  the  outer  discharge  radius  rather  than  a  computed 
average  radius. 
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Additional ly,  this  project  had  as  an  objective  a  significant  expan- 
sion of  the  operating  ranges  examined  beyond  those  previously  considered. 

Because  the  evaluation  of  turbine  performance  parameters  involves  a 
formidable  computation  requirements,  the  data  obtained  in  the  experiments 
were  reduced  by  means  of  the  computer  programs  described  herein,  which  are 
written  in  FORTRAN  IV  language  for  the  IBM  360  computer  system  used  at 
the  Naval  Postgraduate  School. 

The  author  gratefully  acknowledges  the  guidance  and  counsel  of 
Professor  Robert  D.  Zucker,  of  the  Department  of  Aeronautics,  Naval 
Postgraduate  School,  and  of  Professor  Michael  H.  Vavra  of  the  same 
department.   Thanks  are  due  also  to  Mr.  James  E.  Hammer  and  his  coworkers 
on  the  technical  staff  of  the  Turbopropulsion  Laboratory  for  their  many 
hours  of  hard  work  in  modifying,  repairing,  and  operating  the  turbine 
test  rig. 
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SECTION  2 

THE  l.C. P.  RADIAL  TURBINE 

The  turbine  used  in  this  investigation  is  a  dual  discharge,  radial 
inflow  turbine  installed  at  the  Turbopropulsion  Laboratory  of  the  Depart- 
ment of  Aeronautics 9  Naval  Postgraduate  School .   It  was  originally  de- 
signed and  built  by  the  Aerojet  General  Corporation,  for  application  as 
the  power  supply  for  a  missile  propulsion  combustion  cycle.   The  origi- 
nal performance  evaluation  was  conducted  by  Vavra  (1),  with  subsequent 
investigations  by  Finn  (2),  Riley  (3) s  and  Boshoven  (4). 

2.1  Description  of  Turbine  Assembly 

The  basic  turbine  test  installations  used  are  shown  in  Figures  1 
and  2c   The  turbine  assembly  consists  of;   the  wooden  inlet  casing,  with 
its  inner  contour  in  the  shape  of  a  varying  diameter  torus,  or  scroll;  a 
circular  stainless  steel  guide  vane  assembly  consisting  of  seven  circular- 
arc  blades  and  the  rings  themselves;  the  rotor  assembly;  and  two  plexi- 
glass shrouds. 

The  inner  contour  of  the  inlet  casing  was  formed  by  casing  plaster 
of  paris  around  a  wooden  insert  of  the  desired  scroll  dimensions.   The 
inner  surface  of  the  scroll  was  varnished  to  prevent  erosion.   The  wooden 
insert  used  in  this  process  is  shown  in  Figure  3,  and  Figure  4  presents 
a  cross-section  of  the  wooden  casing,  showing  the  shape  of  the  scroll  and 
the  locations  of  the  various  scroll  pressure  taps.   It  is  worthy  of  note 
that  the  flow  entry  is  not  tangential,  as  is  customarily  the  case  in 
radial  turbine  installations.   This  characteristic  is  a  result  of  design 
requirements  in  the  originally  proposed  missile  installation. 
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A  cross-section  of  the  turbine  is  shown  in  Figure  5.   In  this  view 
are  indicated  the  turbine  dimensions  and  minimum  clearance  as  they 
existed  for  the  test  runs.   The  clearance,  or  axial  distance  between  the 
rotor  and  the  shrouds ,  was  increased  during  the  tests  by  means  of  circu- 
lar metal  shims  inserted  between  the  shroud  flanges  and  the  aluminum 
rings  on  which  they  rest. 

The  series  of  tests  conducted  utilized  two  different  rotors.   The 
"dummy"  rotor,  so  called  because  it  does  not  rotate,  is  a  special  device 
designed  to  permit  measurement  of  the  torque  generated  by  the  incoming 
flow  and  thus  provides  a  measure  of  scroll  and  guide  vane  performance. 
This  rotor  is  shown  in  Figure  6.   It  has  a  diameter  of  9.50  inches,  an 
axial  length  of  8.50  inches,  and  36  meridional  blades.   This  diameter 
places  the  blade  leading  edges  at  the  radius  which  corresponds  to  that  of 
the  annulus  containing  the  static  pressure  taps  ahead  of  the  rotor,  and 
the  axial  length  extends  the  blade  trailing  edges  beyond  the  shrouds. 

The  "active"  rotor,  or  the  power-producing  rotating  element,  is 
shown  in  Figure  7,   It  is  9.40  inches  in  diameter  and  has  an  axial  length 
of  5.00  inches.   The  rotor  discharge  radius  is  1.76  inches  at  the  hub  and 
2.94  inches  at  the  tip.   The  fifteen  blades  have  a  spacing  at  the  dis- 
charge varying  from  0.737  to  1.232  inches,  hub  to  tip  respectively.   The 
actual  rotor  discharge  area  on  one  side  of  the  rotor  is  6.43  square  inches, 
including  an  effective  flow  area  through  the  radial  gap  between  the  blades 
and  the  shroud  taken  at  the  discharge  plane  (1).  A  view  of  the  left 
rotor  discharge,  showing  the  plexiglass  shroud  and  several  of  the  eight 
shroud  pressure  taps  which  measure  static  pressure  at  the  rotor  inlet, 
is  given  in  Figure  8. 
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2.2  Modifications  to  Turbine  Assembly 

When  the  turbine  assembly  was  dismantled  in  preparation  for  this 
series  of  tests ,  it  was  discovered  that  one  of  the  plexiglass  shrouds 
was  deeply  gouged.   Since  removal  of  enough  material  to  eliminate  this 
fault  would  have  destroyed  the  contour  of  the  shroud,  it  was  decided  to 
cut  the  contours  of  both  shrouds  so"that  they  would  exactly  fit  the  rotor 
contour.  Additionally,  it  was  decided  to  correct  several  dimensional 
discrepancies  noted  in  earlier  reports. 

Riley  (3)  indicates  a  distance  of  0.943  inch  between  the  inner  ex- 
tremities of  the  shrouds.   Since  the  axial  width  of  the  blade  tips  is 
0.894  inch,  this  indicates  a  basic  (unshimmed)  axial  tip  clearance  of 
0.025  inch  on  each  side  of  the  blade.   Boshoven  (4),  however,  determined 
the  distance  between  shrouds  to  be  0.966  inch,  with  a  resulting  unshimmed 
clearance  of  0.036  inch.  Additionally,  Boshoven  discovered  that  the  in- 
side extremity  of  the  left  shroud  (as  viewed  from  the  turbine  inlet)  was 
offset  0.016  inch  more  from  the  guide  vane  ring  than  the  right  shroud. 

In  order  to  correct  these  discrepancies,  an  extensive  series  of 
measurements  was  undertaken.   The  guide  vane  rings,  which  were  found  to 
be  warped,  were  lapped  until  they  were  flat,  and  all  rough  edges,  burrs, 
and  high  spots  on  the  inner  surfaces  of  the  casing  halves  were  carefully 
removed.   This  resulted  in  a  consistency  to  within  0.0025  inch  in  the 
dimension  between  the  shroud  support  rings,  at  all  points  around  the  pe= 
riphery.   Initial  measurements  showed  a  maximum  variation  of  0.01  inch 
in  this  dimension. 

Measurements  with  the  reworked  shrouds  installed  then  showed  a 
minimum  clearance  of  0.0477  inch,  with  the  left  shroud  offset  axially 
outward  0.0056  inch  more  than  the  right.   The  shroud  flange  thicknesses 
were  then  reduced  so  as  to  compensate  for  this  dissymmetry  and  to  yield 
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a  final  unshimmed  clearance  of  0.015  inch.   Subsequent  measurement  of  the 
axial  travel  of  the  active  rotor9  from  contact  at  one  shroud  face  to  con- 
tact at  the  otherB  verified  that  all  calculated  clearances  used  in  these 
tests  were  correct  to  within  +  0.002  inch. 

Additional  modifications  were  made  to  the  dummy  rotor  in  order  to 
correct  the  axial  stagger  and  unequal  widths  of  the  blades  and  to  correct 
as  much  as  possible  the  large  gaps  which  were  found  to  exist  between  the 
rotor  blade  contours  and  the  shroud  contours.  By  careful  machining,  the 
blade  contours  were  made  uniform9  and  a  much  closer  fit  was  made  to  the 
shroud  contours ,  even  though  it  was  not  possible  to  completely  correct 
this  condition  and  still  maintain  the  blade  width. 

It  is  thus  apparent  that  the  results  obtained  from  the  turbine  after 
this  series  of  modifications  could  differ  to  a  marked  degree  from  those 
previously  obtained  under  identical  operating  conditions.   It  also  seems 
to  this  author  that  both  Riley  (3)  and  Boshoven  (4)  were  operating  at 
clearances  somewhat  greater  than  stated  in  their  respective  studies. 
Based  only  on  these  modifications,  the  distance  between  the  shroud  inner 
extremities  for  their  tests  would  have  been  approximately  0.989  inch, 
yielding  a  minimum  clearance  of  0.048  inch.   This  may  be  compared  with 
the  minimum  clearance  of  0.0245  inch  reported  by  Riley  and  0.036  inch 
reported  by  Boshoven. 

Other  modifications  to  the  components  of  the  test  rig  were  made  and 
are  covered  in  the  sections  to  which  they  apply.   In  this  section  only 
those  modifications  basic  to  the  radial  turbine  itself  are  discussed. 
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ANALYSIS  OF  SCROLL  AND  GUIDE  VANE  PERFORMANCE 

In  order  to  evaluate  the  performance  parameters  of  a  turbine,  it 
is  necessary  to  divide  the  losses  into  those  which  occur  in  the  inlet 
guide  vanes  and  those  which  occur  in  the  rotor  itself.  Riley  (3)  and 
Boshoven  (4)  have  used  the  special  dummy  rotor  installation  previously 
described  (Figure  2)  to  determine  the  losses  in  the  stator,  or  scroll 
and  guide  vane  assembly.  The  installation  determines  the  average  condi- 
tions at  the  rotor  inlet  and  thus  provides  an  indication  of  scroll  and 
guide  vane  efficiency  as  well  as  the  absolute  velocity  (in  both  mangi- 
tude  and  direction)  at  the  rotor  inlet. 

Riley  (3)  concluded  that  the  velocity  coefficient,  which  is  a  measure 
of  the  losses  occurring  in  the  inlet  manifold  and  guide  vanes,  as  well  as 
the  absolute  inlet  flow  angle,  was  essentially  constant  in  the  range  of 
turbine  pressure  ratios  and  axial  clearances  examined.   The  results  ob- 
tained by  Boshoven  (4)  indicate  that  this  assumption  is  invalid;  in 
fact,  the  velocity  coefficients  were  found  to  increase  with  pressure  ratio 
for  a  given  clearance,  and  the  absolute  rotor  inlet  flow  angles  were  found 
to  decrease  linearly  with  pressure  ratio.   Specific  results  of  these  ear- 
lier tests  are  discussed  in  detail  in  Section  3.3. 

It  was  originally  intended  to  utilize  the  dummy  rotor  test  data  of 
Boshoven  (4)  in  this  experiment.  However,  because  of  the  extensive  modi- 
fications previously  discussed,  it  was  necessary  to  repeat  these  tests. 
This  repetition  yielded  the  additional  advantage  of  permitting  testing 
over  a  wider  range  of  clearances. 
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3.1  Instrumentation  and  Test  Procedures 

The  dummy  rotor  shown  In  Figure  6  is  mounted  on  self-aligning  ball 
bearings  and  replaces  the  turbine  rotor.  Its  purpose  is  to  turn  the  flow 
from  the  guide  vanes  into  the  axial  direction  and  to  discharge  it  with  no 
peripheral  velocity  component.  If  this  condition  Is  met,  the  losses  in 
the  scroll  and  guide  vanes  may  be  determined  by  measuring  the  torque  ex- 
erted on  the  static  dummy  rotor,  the  flow  rate,  the  conditions  at  the 
manifold  inlet,  and  the  static  pressure  at  the  dummy  rotor  inlet. 

As  was  previously  mentioned,  the  dummy  rotor  has  an  axial  length  of 
8.50  inches,  which  extends  the  blade  trailing  edges  beyond  the  shrouds. 
To  insure  that  the  discharge  flow  had  no  whirl  components,  flow  straight- 
eners,  also  shown  in  Figure  6,  were  used.  These  devices  consist  of  alumi- 
num caps  containing  a  1.00  inch  depth  of  0.125  by  0.0015  inch  honeycomb 
material.  The  caps  slide  over  the  rotor  shaft  and  cover  the  rotor  blade 
trailing  edges.  The  entire  dummy  rotor  assembly  is  supported  by  brackets 
attached  to  the  manifold  casing,  as  shown  by  Figure  2. 

Air  for  turbine  operation  was  supplied  by  an  Allis-Chalmers  12-stage 
axial  compressor.  Air  from  the  main  supply  line  passes  through  a  four- 
Inch  pipe,  a  settling  tank,  a  five-inch  pipe,  and  then  into  the  turbine. 
The  flow  rate  was  regulated  by  remotely  controlled  butterfly  valves  in  the 
inlet  line  to  the  turbine  and  in  the  outside  discharge  line  of  the  Allis- 
Chalmers  compressor.   (See  Figure  9.)  The  controls  for  operating  these 
valves  are  located  in  the  control  room  (Figure  10). 

The  flow  rate  of  the  incoming  air  was  measured  with  a  sharp-edged 
orifice  located  in  the  four-inch  inlet  pipe.  This  orifice  is  also  shown 
in  Figure  9.  The  pressure  ahead  of  the  orifice  and  the  pressure  differ- 
ential across  the  orifice  were  measured  with  standard  vena  e en tracts 
taps.  The  total  temperature  ahead  of  the  orifice  was  measured  with  an 
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iron-constantan  thermocouple  installed  in  a  Kiel  temperature  probe.   The 
orifice  installation  conforms  to  standards  set  forth  by  Stearns,  e_t  al. 
(5). 

The  total  inlet  conditions  in  the  five-inch  pipe  just  upstream  of 
the  turbine  were  measured  with  a  Kiel  probe „  also  incorporating  an  iron- 
constantan  thermocouple.   The  static  pressure  at  the  same  location  was 
measured  by  two  static  pressure  taps  in  the  line.   The  static  pressure 
at  the  rotor  inlet  was  measured  by  pressure  taps  located  at  eight  pe- 
ripheral stations  in  each  shroud.   The  axial  location  of  these  taps  is 
depicted  in  Figure  59  and  their  peripheral  locations  are  shown  in  Figure 
4.   These  taps  were  moved  to  the  positions  shown  from  their  original  loca- 
tions near  the  shroud  inner  extremities  by  Boshoven  (4)  in  order  to  re- 
duce the  effects  of  local  pressure  perturbations  experienced  in  previous 
tests. 

The  torque  exerted  on  the  dummy  rotor  by  the  flow  was  measured  on  a 
25-pound  capacity  Toledo  precision  scale 9  shown  in  Figure  2.   A  twelve- 
inch  lever  arm  was  attached  to  the  rotor  shaft  and  incorporated  a  simple 
provision  for  leveling9  so  that  the  force  was  exerted  on  the  scale  through 
an  exactly  vertical  link.   This  was  done  to  preserve  the  direct  relation- 
ship between  scale  reading  and  moment.  A  counter  balance  was  added  to  the 
lever  arm  in  order  to  extend  the  capacity  of  the  scales  to  accommodate 
higher  mass  flow  rates. 

A  general  schematic  of  the  instrumentation  used  in  this  and  the  sub- 
sequent active  rotor  tests  is  shown  in  Figure  11.   From  this  schematic  it 
may  be  seen  that  the  pressure  upstream  of  the  flow  measuring  orifice, 
the  total  and  static  pressures  at  the  turbine  inlet,  and  the  reference 
pressures  1  and  2  were  read  from  an  eleven  column,  40-inch  mercury  mano- 
meter board  located  in  the  control  room.   This  board  was  referenced  to 
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atmospheric  pressure.  Additionally,,  one  arbitrarily  selected  scroll  pres- 
sure tap  was  connected  to  the  mercury  board,  and  this  pressure,  together 
with  the  total  inlet  pressure,  was  used  to  set  the  pressure  ratio  across 
the  scroll  and  guide  vanes.   The  total  inlet  pressure  displayed  on  the 
mercury  board  was  used  only  for  setting  this  pressure  ratio  since  a  more 
representative  total  pressure  was  computed  in  the  data  reduction. 

The  sixteen  shroud  pressure  taps  that  measure  the  static  pressure  at 
the  rotor  inlet  were  read  on  a  fifty-column,  96-inch  water  manometer  board. 
Also  displayed  on  this  board  were  the  pressures  at  various  stations  in  the 
scroll  itself,  although  these  pressures  were  not  used  in  the  tests.   These 
pressures  were  measured  against  atmospheric  pressure  for  the  entire  series 
of  dummy  rotor  runs.   However,  additional,  controlled,  pressures  could  have 
been  imposed  on  the  water  board  reservoirs  by  means  of  regulated  house  air 
if  this  had  become  necessary  in  order  to  avoid  exceeding  the  limits  of  the 
manometer  board.   These  pressures  (the  reference  pressures  1  and  2  men- 
tioned above)  were  displayed  on  the  mercury  manometer  board.   The  first  of 
these  pressures,  displayed  on  column  10  of  the  mercury  board,  was  recorded 
as  a  test  datum  since  it  applied  to  the  sixteen  columns  which  indicated 
the  rotor  inlet  pressures.   The  value  of  this  pressure  was  thus  zero  for 
all  dummy  rotor  runs. 

The  pressure  differential  across  the  flow  orifice  was  measured  on  a  20- 
inch  Meriam  model  34FB2  mercury  micro-manometer,  which  has  a  measuring  ac- 
curacy of  +0.001  inch.   A  pair  of  these  instruments  is  shown  in  Figure  12. 

The  potentials  generated  by  the  iron-constantan  thermocouples  were 
measured  on  a  48-channel  Brown  potentiometer,  manufactured  by  Minneapolis- 
Honeywell.   This  potentiometer,  shown  in  Figure  13,  was  located  in  the 
control  room  and  used  an  ice  bath  in  the  test  cell  as  a  cold  junction 
reference. 
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A  series  of  eighteen  test  runs  were  made  on  the  dummy  rotor  instal- 
lation, of  which  eight  are  documented  in  this  report.   These  eight  runs 
represent  clearances  of  0.015,  0.024,  0.030,,  0.041,  0,051,  0.061,  0.072, 
and  0.081  inches ,  respectively.   For  each  clearance  ten  pressure  ratios 
were  sampled,  ranging  in  value  from  1.05  to  1.65.   These  pressure  ratios 
represent  the  ratio  of  total  turbine  inlet  pressure  to  the  static  pres- 
sure at  the  rotor  inlet. 

This  large  number  of  runs  was  necessitated  by  several  factors.   Data 
from  the  initial  runs  possessed  an  unacceptable  degree  of  scatter.   This 
was  attributed  to  an  unstable  supply  of  source  air  from  the  Allis-Chalmers 
compressor.   When  this  unsatisfactory  trend  was  discovered,  the  speed 
control  system  was  overhauled,  and  the  subsequent  test  results  were  signi- 
ficantly improved. 

Additionally,  the  stability  of  the  system  was  improved  by  judicious 
use  of  the  aforementioned  butterfly  control  valve.   Experimentation  showed 
that  if  less  of  the  excess  air  from  the  compressor  were  discharged  to 
the  outside,  the  compressor  would  be  operating  at  a  higher  pressure  ratio, 
(with  a  greater  pressure  drop  across  the  valve  feeding  the  turbine) ,and 
the  pressurized  air  in  the  house  plumbing  would  act  as  an  additional  plenum. 
The  system  fluctuations  would  then  be  substantially  reduced.   These  im- 
provements in  equipment  and  operating  technique,  as  well  as  the  inevita- 
ble improvement  in  data  taking  technique,  led  to  such  improved  results  that 
many  runs  were  repeated  in  order  to  obtain  the  best  possible  test  data. 

As  first  observed  by  Boshoven  (4),  it  was  noted  that  the  pressure 
distribution  around  the  shroud  was  very  sensitive  to  the  relative  posi- 
tions of  the  rotor  blade  leading  edges.   To  eliminate  this  factor  as  a 
possible  variable,  it  was  arranged  to  permit  the  return  of  the  rotor  to 


the  same  location  each  time  a  new  pressure  ratio  was  set.   This  was  ac- 
complished by  means  of  a  wing  nut  on  the  end  of  the  rod  connecting  the 
moment  arm  to  the  scale. 

It  was  also  observed  that  the  pressures  seemed  to  be  higher  by  about 
4  to  5  inches  of  water  at  the  right  side  of  the  rotor  inlet  than  at  the 
left.   This  condition  could  be  eliminated  by  moving  the  rotor  approximately 
0.001  inch  axially  to  the  left;  but  since  the  condition  was  also  present  in 
the  scroll  itself,  it  was  attributed  to  scroll  pressure  distribution  rather 
than  to  misalignment.   The  rotor  was  centered  for  all  runs. 

3.2  Theoretical  Analysis 

From  the  law  of  moment  of  momentum  (6),  the  moment  M  exerted  on  the 
dummy  rotor  by  the  flow  is  given  by 

(v\-(W/g)R1Vu1  a) 

where  W  is  the  mass  flow  rate,  R  is  the  inlet  radius  of  the  dummy  rotor, 
and  V   is  the  peripheral  velocity  component  at  R. .   Equation  (1)  is 
valid  for  the  case  of  steady  flow  in  the  dummy  rotor  because  the  pressures 
acting  over  the  inlet  and  discharge  areas  cannot  produce  a  moment  about 
the  axis,  and  because  the  flow  straighteners  insure  an  axial  discharge 
velocity. 

A  measure  of  the  performance  of  the  scroll  and  guide  vanes  is  given 
in  the  form  of  the  velocity  coefficient,  (p   ,  which  is  defined  by 

where  V  is  the  average  velocity  at  the  dummy  rotor  inlet  and  V    is 
the  theoretical  velocity  at  the  dummy  rotor  inlet  which  would  result  from 
an  isentropic  expansion  from  the  manifold  inlet  total  conditions  P.~>  T  ~, 
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to  the  static  pressure  p  at  the  dummy  rotor  inlet.   This  process  is 
shown  on  the  T-s  diagram  of  Figure  14.   The  component  velocities  may  be 
seen  on  Figure  15. 

From  the  energy  equation  for  an  adiabatic,  no  work  process,  the 
velocity  V1  at  the  dummy  rotor  inlet  is 

Vi  -  /£gTcp  (Tto-Tt)'  O) 

If  use  is  made  of  the  isentropic  pres sure - temperature  relation  for  a 

perfect  gas  and  of  Equation  (2) ,  V.  may  be  alternately  expressed  as 


Vi  =  </>  /^gjcPTt0[i-(-pi/p^¥] 


(4) 


Application  of  the  equation  of  continuity  at  the  dummy  rotor  inlet 
yields  the  following  expression  for  the  average  meridional  velocity 
component: 

Vmi  =  W/  fi  Ai  (5) 

where  a>  is  the  density  at  the  dummy  rotor  inlet,  and  A  is  the  cross- 
sectional  area  normal  to  the  flow  (meridional  cross-section).  The  area 
A.  is  given  by 

Ai  =    Zff  RL  Bi  <«> 

where  B,  is  the  distance  between  the  shrouds  at  the  rotor  inlet.   Using 

i 

the  equation  of  state  for  a  perfect  gas,  £>  may  be  expressed  as 

ft  =  Tt/RgTi  m 

Substitution  of  this  relation  and  Equation  (6)  into  Equation  (5)  yields 
the  resultant  expression  for  ¥  . : 
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From  Equation  (1) ,  the  peripheral  velocity  may  be  expressed  as 


and  V-  may  additionally  be  given  by 


Vt  -  /Vm*  +  V-vu*' 


(10) 


Thus,  three  equations  expressing  the  velocity  V.  are  established, 
these  being  Equations  (3),  (4),  and  (10).   Upon  substitution  of  Equations 
(8)  and  (9)  into  Equation  (10),  it  may  be  seen  that  this  set  of  three 
equations  contains  three  unknowns.   The  unknown  quantities  to  be  deter- 
mined are  the  velocity  V-  ,  the  velocity  coefficient  (p  ,  and  the  tempera- 
ture T  ,  all  other  quantities  being  either  constants  or  determined  direct- 
ly from  measured  data.   Solution  of  this  system  of  equations  establishes 
the  velocity  components,  which  may  then  be  used  to  calculate  the  absolute 
flow  angle  at  the  dummy  rotor  inlet,  according  to  the  relation 

A  nozzle  loss  coefficient,  Ip  ,  may  be  obtained  from  the  relation 

where  t>  is  an  efficiency  based  on  the  ratio  of  the  actual-to-isentropic 
temperature  drop  across  the  scroll  and  guide  vanes.   It  is  introduced  here 
only  to  illustrate  the  significance  of  \o      and  (h . 
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3.3  Reduction  of  Data:   Program  SCROLL 

The  computer  program  SCROLL,  listed  in  Appendix  A,  determines  the 
unknown  quantities,  V. ,  (f>   ,  and  T  ,  and  the  absolute  inlet  angle,  ©<.  , 
from  the  test  data  using  an  iterative  procedure.   The  program  makes 
temperature  corrections  to  the  physical  properties  of  air  and  to  the 
specific  gravities  of  the  water  and  mercury  used  as  manometer  fluids. 
It  converts  the  manometer  readings  and  thermocouple  voltages  into  pres- 
sures and  temperatures ,  as  required s  and  then  establishes  the  mass  flow 
rate.  The  average  conditions  at  the  dummy  rotor  inlet  are  then  calcu- 
lated from  the  measured  torque  in  accordance  with  the  previously  given 
relations. 

Program  SCROLL  has  been  described  in  detail  by  Riley  (3)  and  Boshoven 
(4).   Except  for  conversion  to  the  FORTRAN  IV  language,  the  program  used 
in  this  experiment  remains  essentially  the  same;  therefore,  the  details 
of  the  computational  processes  will  not  be  repeated.   However,  certain 
minor  changes  were  incorporated  in  the  temperature  corrections  applied  to 
the  manometer  fluids,  in  the  thermocouple  equations  used,  and  in  the  pres- 
sure and  flow  rate  solutions.   These  changes  are  briefly  described  in 
Appendix  A,  although  they  become  evident  if  the  program  used  in  this  re- 
port is  compared  with  either  of  the  two  above-listed  references. 

3.4  Discussion  of  Results 

The  results  of  program  SCROLL  consist  primarily  of  values  of  the 
velocity  coefficient,  d   ,  and  the  absolute  flow  angle  at  the  rotor  in- 
let, ©<  •   These  values  are  plotted  versus  M. ,  the  Mach  number  at  the 
rotor  inlet,  as  opposed  to  the  pressure  ratio,  P  /p  ,  because  the  Mach 
number  is  a  more  significant  parameter  in  the  application  of  test  results 
to  geometrically  similar  machines.   Figure  16  shows  <p     as  a  function  of 
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M  for  each  of  the  eight  clearances  tested,  and  Figure  17  is  a  similar 
plot  for  o<  as  a  function  of  M.. .   Each  of  the  curves  was  fitted  to  the 
data  obtained  for  that  run  by  the  least  squares  method.   The  tabulated 
results  for  each  of  the  eight  runs  are  given  in  Appendix  A. 

It  may  be  seen  from  Figure  16  that  the  velocity  coefficient 
increases  with  Mach  number,  the  average  increase  being  approximately  4.75 
per  cent  over  a  Mach  number  range  of  0.2  to  about  0.8.  This  percentage 
increase  is  practically  the  same  regardless  of  clearance.   This  increase 
is  due,  at  least  in  part,  to  the  decreasing  coefficient  of  friction  in 
the  scroll  and  guide  vanes,  which  results  from  the  increasing  Reynolds 
number. 

The  coefficient  (j>   generally  increases  with  decreasing  clearance  for 
a  given  pressure  ratio  with  the  maximum  values  of  d)   occurring  at  the  mini- 
mum clearance  of  0.015  inches.  As  clearance  is  increased  beyond  a  value 
of  about  0.070  inches,  the  incremental  drop  in  (f>    per  increment  of  clear- 
ance increase  appears  to  become  much  larger.  An  explanation  for  this  oc- 
currence is  that  as  clearance  is  increased,  an  increasing  amount  of  flow 
escapes  around  the  blade  tips.   This  flow  loses  angular  momentum  by  wall 
friction  and  subsequently  strikes  the  dummy  rotor  blades  at  some  lesser 
radius,  thus  imparting  a  reduced  amount  of  torque  to  the  rotor.   The  cumu- 
lative effect  of  this  energy  loss  would  logically  become  more  pronounced 
at  increased  clearances.   The  overlapping  of  certain  of  the  curves  in  Fig- 
ure 16  at  the  lower  Mach  numbers  occurred  because  of  the  relative  in- 
sensitivity  of  the  instrumentation  at  low  pressures. 

Figure  17  shows  that  the  absolute  rotor  inlet  flow  angle,  oC,,  in- 
creases with  increasing  Mach  number  up  to  a  Mach  number  of  approximately 
0.4.   Beyond  that  point,  oC.  decreases  with  increasing  Mach  number.   These 
changes  are  small,  although  they  appear  to  be  significant  because  of  the 
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highly  expanded  scale  of  Figure  17.   In  fact,  the  average  maximum  change 

in  o<.  over  the  range  of  Mach  numbers  is  approximately  0.52  degrees.   In 

general ,  increasing  the  mass  flow  rate  will  increase  V1  and,  of  course, 

the  components  V  -  and  V  , .   The  variations  in  the  magnitudes  of  these 
ul      ml  ° 

two  quantities  as  the  Mach  number  is  increased  determine  the  variations 
of  c<  ,  as  shown  by  Equation  (11). 

It  is  also  shown  on  Figure  17  that  the  angle  cK.  increases  as  clear- 
ance increases.   This  is  primarily  because  V  ,  is  inversely  proportional 

ml 

to  the  meridional  cross  section,  which  increases  with  clearance  since  it 

depends  on  the  distance  between  the  shrouds.   By  Equation  (11),  the  smaller 

values  of  V  ,  result  in  larger  inlet  flow  angles, 
mi 

It  was  concluded  by  Riley  (3)  that  the  velocity  coefficient  <j>   was  es- 
sentially independent  of  the  pressure  ratio  P  /p  ,  and  that  a  constant 
representative  value  of  0.889  was  sufficiently  accurate.  Although  he 
determined  that  &(.    decreases  with  increasing  pressure  ration,  he  also 
assumed  a  representative  value  of  80.0  degrees  for  o<-,  • 

Boshoven  (4)  subsequently  investigated  the  scroll  and  guide  vane  per- 
formance in  greater  detail  and  obtained  results  that  agree  in  a  qualita- 
tive sense  with  those  described  above.   Because  of  the  many  modifications 
made  to  the  machine  components,  both  in  preparation  for  this  experiment 
and  during  the  course  of  Boshoven "s,  a  more  precise  comparison  is  not 
warranted.   It  is  worthy  of  note,  however,  that  Boshoven  did  not  obtain 
the  continuous  decrease  in  <fi   with  increasing  clearance  which  is  apparent 
from  Figure  16.   Further,  the  increase  in  <j>   with  Mach  number  found  in  the 
present  tests  is  more  pronounced  than  Boshoven's,  his  increase  being  on 
the  order  of  2   per  cent  for  a  given  clearance,  over  the  range  of  Mach 
numbers.   It  is  felt  that  these  changes  have  occurred  because  of  the 
machining  of  the  rotor  so  that  it  more  closely  fits  the  shroud  contours, 
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thus  minimizing  flow  loss  around  the  blades  and  permitting  closer  mini- 
mum clearances. 

An  examination  of  the  data  points  given  in  Figure  16  for  0  as  a 
function  of  Mach  number  shows  that  some  scatter  exists.   This  scatter 
may  be  attributed  to  fluctuations  in  the  air  supply  from  the  Allis- 
Chalmers  compressor.   The  presence  of  this  problem  and  methods  under- 
taken to  alleviate  it  have  been  discussed  previously  in  Section  3.1. 
This  scatter  of  data  is  also  present  in  the  curves  of  c<  as  a  function 
of  Mach  number  given  in  Figure  17,  but  to  a  lesser  extent. 

In  order  to  use  the  results  of  the  dummy  rotor  tests  in  later  data 
reduction  programs,  it  was  necessary  to  obtain  equations  representing  0 
and  oC   in  the  form 

©C   =  Constant  [f   (clearance) ] [f   (Mach  number)] 
The  function  f.  was  obtained  by  first  taking  the  ratio  of  oC     to  of 
maximum  at  each  Mach  number,  for  each  clearance.   The  values  of  this 
ratio  were  then  plotted  versus  clearance.   From  this  plot  an  average  curve 
was  obtained,  and  its  equation  was  computed  by  the  least  squares  technique, 
The  resulting  equation  expressed  °^,  /  °C,  as  a  dimensionless  function 

of  clearance. 

The  function  f _  was  obtained  by  computing  the  ratio  of  ©^  to  the 
value  of  tf^C.  at  a  Mach  number  of  0.5,  and  plotting  this  ratio  versus 
Mach  number  for  each  clearance.   Again,  an  average  curve  was  determined, 
and  its  equation  was  computed.   This  equation,  again  dimensionless,  ex- 
presses o£- 1  o(.  ,M_  ,-\  as  a  function  of  Mach  number. 

The  appropriate  scaling  factor  was  determined  for  each  test  point. 
An  average  of  all  of  these  factors,  which  did  not  differ  greatly  in  value, 
was  taken  as  the  constant  in  the  above  expression.   The  equation  for  0 
was  determined  in  the  same  manner. 
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The  final  equations  resulting  from  this  procedure  were 

°<i.  =  {&.l(ol)l.%ZC,l   -I-  .2W8  (Cine)  "  •  &50VO   (Cine)  Z   ] 

(13) 


<fi=   (•  92173) [i. 01^33-  LZ4fc0i  (Cine)  4  24.0£8?9  (CUc)x 

-  \S0.ZCZ5Q  (Cine)11  J 
[  .  93i49  f  .  !8Z9z(«I)-  .O^QZC/Hj)*] 


Families  of  curves  produced  by  these  equations  are  given  in  Figures  18 
and  19,  for  <p   and  <X-    respectively.   These  curves  agree  well  with  those 
of  Figures  16  and  17  when  consideration  is  given  to  the  approximations 
introduced  in  the  derivations  of  Equations  (13)  and  (14).  For  the  veloc- 
ity coefficient,  ^  ,  the  maximum  difference  which  exists  at  any  point 
between  the  curves  of  Figure  16  and  those  of  Figure  18  is  0.0052.  A 
similar  comparison  for  o(     yields  a  maximum  difference  of  0.06  degrees. 
Because  of  the  manner  in  which  Equations  (13)  and  (14)  were  derived,  the 
most  significant  discrepancies  occur  near  the  extremes  of  the  Mach  number 
range . 
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SECTION  4 

ANALYSIS  OF  OVERALL  TURBINE  PERFORMANCE 

The  turbine  overall  performance  tests,  when  combined  with  the  results 
of  the  dummy  rotor  tests,  provided  data  from  which  detailed  performance 
parameters  of  the  turbine  could  be  evaluated.   The  specific  parameters 
evaluated  were:   the  total-to-static  overall  efficiency;  velocity  and 
loss  coefficients  for  the  rotor  and  guide  vanes;  rotor  inlet  and  outlet 
flow  angles;  referred  speed,  moments,  and  flow  rates;  the  head  coefficient; 
and  the  degree  of  reaction.   Calculations  were  made  assuming  no  bearing 
losses  and  also  for  the  case  of  "maximum"  bearing  losses.   Bearing  loss 
values  were  determined  from  equations  given  by  Vavra  (1),  which  were 
established  from  coast-down  tests.   The  analysis  was  made  using  a  "mean" 
streamline  approach,  which  is  defined  as  utilizing  a  representative 
streamline  along  which  the  measured  flow  conditions  are  assumed  to  exist. 
This  streamline  was  taken  as  that  at  the  outer  discharge  radius  of  the 
turbine  for  this  analysis,  and  thus  all  losses  are  referred  to  that  outer 
radius. 

4.1   Instrumentation 

The  "active"  rotor  shown  in  Figure  7  was  installed  for  this  series 
of  tests,  and  the  rotor  and  wooden  casing  were  placed  in  an  aluminum 
test  stand.  Figure  1  shows  the  test  installation.  Alignment  of  the  scroll 
with  respect  to  the  rotor  was  accomplished  by  the  use  of  the  dummy  align- 
ment shaft  and  dial  micrometers  shown  in  Figure  20.   Precise  alignment 
was  found  to  be  quite  difficult  because  the  massive  wooden  casing  and 
scroll  assembly  had  to  be  positioned  relative  to  the  rotor.   Positioning 
was  accomplished  by  means  of  four  adjusting  bolts  located  at  the  bottom 
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of  the  casing  and  by  four  angle-irons  attached  to  the  casing  and  resting 
on  the  bearing  supports.   In  spite  of  the  difficulty  involved,  alignment 
of  the  concentric  aluminum  rings  which  serve  as  supports  for  the  shrouds 
was  obtained  to  within  +  0.001  inch  in  the  radial  direction  from  the 
dummy  alignment  shaft.   These  rings  were  found  to  deviate  from  a  plane 
perpendicular  to  the  rotor  axis  by  less  than  0.002  inch  at  a  radius  of 
4.7  inches. 

The  power  developed  by  the  turbine  was  absorbed  by  a  Vortec  Products 
Company  air  dynamometer9  model  20.075 9  shown  in  Figure  1.   Lubrication  of 
the  dynamometer  was  accomplished  by  an  oil  mist  generated  by  air  pressure. 
The  torque  output  was  measured  with  a  torque  capsule  having  a  maximum 
rated  capacity  of  400  inch-pounds.   This  capsule  utilized  strain  gages 
arranged  on  an  internal  flexure  and  was  additionally  instrumented  with  an 
iron-cons tantan  thermocouple  to  provide  a  means  of  monitoring  the  capsule 
temperature.   The  signal  generated  by  the  strain  gages  was  read  on  a  Day- 
tronic  model  700  digital  indicator  located  in  the  control  room  (Figure  10). 

The  dynamometer  shroud  was  modified  from  the  configuration  used  by 
Riley  (3)  so  as  to  discharge  the  hot  exhaust  air  radially  outward,  rather 
than  back  onto  the  torque  capsule  and  housings,  in  order  to  minimize  ther- 
mal effects  on  the  torque  capsule.   The  circular  aluminum  ring  visible  in 
Figure  1  was  used  to  accomplish  this  modification.   Additionally,  a  wooden 
baffle,  not  shown  on  Figure  19  was  used  to  protect  the  capsule  from  the 
cold  turbine  discharge  air9  and  to  prevent  any  interference  between  the 
dynamometer  discharge  and  the  right  rotor  discharge. 

Figure  21  shows  the  inner  face  of  the  torque  capsule,  the  steel  quill 
shaft  used  to  connect  the  dynamometer  to  the  turbine ,  and  the  magnetic  pick- 
up and  six-lobe  flux  cutter  used  to  measure  turbine  speed.   The  turbine 
speed  was  read  from  a  Hewlett-Packard  electronic  counter  located  in  the 
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control  room  (Figure  22).   Speed  was  regulated  by  means  of  a  remotely 
controlled  motor  on  the  dynamometer  itself,  which  changes  the  load  capac- 
ity of  the  dynamometer.   The  remote  control  unit  is  shown  in  Figure  10, 
and  the  drive  motor  is  shown  in  Figure  1. 

Lubrication  of  the  turbine  bearings  was  provided  by  a  dry-sump  system 
utilizing  a  gear-driven  oil  pump  and  an  upper  and  a  lower  oil  mist  injector 
for  each  bearing.   The  lubricating  fluid  used  was  1010  jet  engine  oil. 
The  cart  containing  the  pump,  oil  reservoir,  heat  exchanger,  and  system 
pressure  gages  is  shown  in  Figure  23,  and  the  lubrication  system  is  shown 
schematically  in  Figure  24.  For  these  tests  the  previously  existing  sys- 
tem was  changed  to  that  shown  in  Figure  24  to  permit  measurement  of  the 
oil  temperature  into  the  bearings  and  at  the  point  of  discharge  from  both 
bearings.   This  was  accomplished  with  iron-constantan  thermocouples  in- 
stalled as  shown  on  the  schematic. 

Figure  24  also  shows  the  iron-constantan  thermocouples  used  to  mea- 
sure the  bearing  outer  race  temperatures.  Vibration  signals  from  the 
bearings  were  detected  by  a  piezo-quartz  accelerometer  mounted  on  the  up- 
per left  bearing  cap  and  monitored  on  a  Panoramic  vibration  analyzer  manu- 
factured by  the  Singer  Metrics  Company.   The  accelerometer  is  visible  in 
Figure  8,  and  the  analyzer  is  shown  in  Figure  22. 

The  pressures  and  flow  rates  were  measured  in  the  same  manner  as  for 
the  dummy  rotor  test  previously  described.   The  signals  from  the  thermo- 
couples were  initially  read  from  the  Brown  potentiometer  shown  in  Figure 
13.   During  the  latter  part  of  the  testing  period,  however,  they  were 
read  from  a  new  automatic  data  acquisition  system  manufactured  by  B  &  F 
Instruments,  Inc.   This  system  is  shown  in  Figure  25. 
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4.2  Test  Procedures 

Test  runs  were  conducted  at  five  different  axial  clearances:  0.081, 

0.061,  0.041,  0.024,  and  0.015  inches.  For  each  clearance,  runs  were 

made  at  total  inlet  to  static  discharge  pressure  ratios  of  1.2,  1.3,  1.4, 

1.5,  1.6,  and  1.7.  At  each  pressure  ratio,  turbine  speed  was  varied  from 

the  minimum  speed  (maximum  load  condition)  for  that  pressure,  to  the 

highest  speed  attainable,  in  increments  of  about  1,500  rpm.   Maximum 

speed  was  restricted  to  18,000  rpm  for  reasons  of  safety. 

Measurement  of  the  torque  at  zero  speed  was  also  made  for  each  pres- 
sure ratio  at  each  clearance.   This  was  accomplished  by  engaging  the  dyna- 
mometer rotor  lock  and  measuring  the  resulting  torque. 

As  turbine  speed  was  increased  during  the  first  run,  a  point  was 
reached  at  which  the  bearings  overheated  and  started  to  seize,  causing 
a  rapid  drop  in  speed.   The  turbine  was  dismantled  and  carefully  examined. 
No  apparent  damage  had  been  incurred  by  the  ball  bearings.   The  oil  sys- 
tem was  flushed  and  refilled,  and  the  unit  was  reassembled.   This  condi- 
tion repeated  itself  in  several  subsequent  runs,  occurring  whenever  the 
bearing  temperature  reached  approximately  135°F.  Additionally,  consider- 
able amounts  of  oil  were  being  lost  from  around  the  shaft  and  bearing 
caps  at  higher  speeds.  Apparently  neither  of  these  conditions  had  ex- 
isted in  previous  tests. 

Since  the  bearings  themselves  did  not  appear  to  be  at  fault,  ex- 
tensive tests  were  made  to  determine  whether  changing  the  amount  of  lubri- 
cating oil  flow  would  alleviate  these  problems.   It  was  found  that  vary- 
ing the  oil  system  pressure  from  10  pounds  per  square  inch  to  50  pounds 
per  square  inch  would  reduce  the  bearing  temperature  a  maximum  of  six 
degrees,  but  this  variation  did  not  significantly  reduce  leakage.  A 
solution  was  ultimately  found  to  both  problems  by  carefully  positioning 
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the  oil  injectors  so  that  they  discharged  directly  into  the  balls  of  the 
bearing  and  did  not  strike  the  race  and  by  slightly  enlarging  the  injector 
orifices  so  as  to  increase  the  mass  flow  rate  of  the  oil.  The  bearing 
temperatures  could  then  be  kept  below  a  maximum  of  about  125 °F  at  18,000 
rpm,  and  oil  leakage  was  negligible. 

The  fact  that  the  bearings  experienced  a  maximum  temperature  in- 
crease of  about  50°  over  the  range  of  speeds  tested,  whereas  in  previous 
tests  the  increase  had  been  only  about  half  as  much  (7),  created  a  ques- 
tion as  to  the  applicability  of  the  bearing  loss  information  obtained 
from  the  coast-down  tests  conducted  by  Vavra  (1).   It  was  therefore  decided 
to  attempt  to  evaluate  the  bearing  losses  by  measuring  the  oil  flow  rate 
and  the  amount  of  heat  absorbed  by  the  oil  as  it  flowed  across  the  bear- 
ings.  If  the  constant  pressure  specific  heat,  C  ,  of  the  oil  were  known, 
the  energy  lost  in  the  bearings  could  be  obtained  from  the  relation 

BFM  -       —  CpoU     A~rb&Ar,^s  (15) 

XjO" 

where  BFM  is  the  bearing  friction  moment  in  foot-pounds  and  ~uj    is  the 
turbine  rotational  speed. 

Numerous  tests  were  made  in  order  to  calibrate  the  mass  flow  rate  of 
the  oil  with  the  system  pressure,  and  the  lubrication  system  was  modified 
to  include  the  previously  discussed  thermocouples  ahead  of  and  after  the 
bearings.   Flow  rate  calibration  was  accomplished  by  weighing  the  amount 
of  oil  discharged  by  the  injectors  on  each  side  in  a  given  time  interval. 

In  attempting  to  evaluate  bearing  losses  by  this  technique,  the  re- 
sults obtained  from  several  test  points  indicated  losses  on  the  order  of 
20  per  cent  of  the  minimum  losses  predicted  by  Vavra  (1).   In  fact,  it 
was  discovered  that  the  exposure  of  the  inlet  lines,  the  bearing  block  it- 
self, and  the  sumps  to  the  cold  discharge  air  sometimes  produced  cooler 
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temperatures  on  the  discharge  side  of  the  bearings  than  existed  ahead 
of  the  injectors .   Since  insulation  of  the  oil  lines  and  the  bearing 
blocks  to  the  extent  necessary  to  overcome  this  deficit  was  not  feasible, 
this  attempt  to  determine  the  bearing  losses  was  abandoned.   The  oil 
temperatures  were  recorded  as  data  for  future  use,  however,  and  appear 
in  Appendix  B. 

The  dynamometer  torque  capsule  was  calibrated  both  before  and  after 
each  day's  runs  in  order  to  permit  consideration  of  any  variations  during 
the  course  of  the  runs.   Calibration  was  accomplished  by  use  of  the  de- 
vice shown  in  Figure  26.   An  aluminum  plug  was  fabricated  to  fit  into  the 
torque  ring  of  the  capsule.   Torque  was  applied  by  loading  the  weight 
pans  on  the  moment  arm,  which  was  attached  to  the  plug.   Riley  (3)  at- 
tached the  moment  arm  directly  to  the  inlet  side  of  the  dynamometer; 
thus,  a  moment  was  introduced  due  to  the  axial  distance  between  the 
capsule  and  the  point  of  application  of  the  load.   The  apparatus  used  in 
these  tests  reduced  this  potential  source  of  error  almost  entirely. 

Loads  of  0  to  400  inch-pounds  were  applied  to  the  capsule  in  100 
inch-pound  increments.   The  span  of  the  digital  readout  was  set  at  0  to 
1000  counts,  the  latter  approximately  corresponding  to  400  inch-pounds. 
Calibration  readings  were  taken  for  both  increasing  and  decreasing  loads 
in  order  to  evaluate  the  amount  of  hysteresis  present  in  the  readout. 
As  each  load  was  applied,  the  capsule  housing  was  tapped  gently  with  a 
mallet  to  overcome  any  drag  in  the  capsule  itself  and  to  eliminate  errors 
caused  by  lag  in  the  servomechanism  of  the  digital  counter.   This  pro- 
cedure is  valid,  because  a  like  amount  of  vibration  was  observed  to  exist 
under  the  dynamic  load  conditions  of  an  actual  test  run.   Using  this 
technique,  the  amount  of  hysteresis  present  between  the  loading  and  un- 
loading cycles  was  found  to  be  negligible. 
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The  torque  readouts  for  calibrations  performed  both  before  and  after 
the  test  runs  were  found  to  vary  linearly  with  applied  torque.   However, 
the  data  obtained  after  the  run  were  generally  higher ,  by  2  to  3  counts 
at  zero  load,  to  13  to  15  counts  at  400  inch-pounds,  than  the  original 
calibration.   This  was  felt  to  be  due  to  the  thermal  effects  of  the  higher 
capsule  temperatures  which  existed  during  the  run.   A  thermocouple  was 
installed  in  the  capsule  in  an  attempt  to  establish  some  correlation  be- 
tween this  temperature  rise  and  the  torque  count  increase.   Tests  in  which 
the  capsule  was  heated  from  the  outside  with  a  strip  heater  were  incon- 
clusive, because  similar  increase  in  torque  counts  for  given  torque  load- 
ings could  not  be  produced,  even  though  the  temperatures  indicated  by  the 
thermocouple  were  considerably  higher  than  those  recorded  during  the  actual 
test  runs. 

When  the  turbine  is  in  operation,  the  capsule  is  heated  from  within 
by  heat  from  the  dynamometer  bearings.   This  is  modified  by  the  dynamometer 
water  cooling  system  and  by  the  cold  turbine  discharge  air.   A  method  of 
reproducing  the  interaction  of  all  of  these  conditions  in  static  calibra- 
tion tests  was  not  found,  and  so  the  torque  calibration  data  ultimately 
used  represented  an  average  of  the  data  obtained  before  and  after  the  test 
runs.   This  accounts  for  some  increase  due  to  thermal  effects  and  is  felt 
to  be  a  reasonable  representation.   In  some  instances,  the  capsule  tempera- 
ture was  observed  to  increase  immediately  following  shut-down,  while  at 
other  times  the  reverse  effect  occurred.   It  is  felt  that  more  representa- 
tive temperature  data  could  be  obtained  by  relocating  the  thermocouple  on 
the  flexure  closer  to  the  strain  gages. 

The  compressor  which  provides  house  air  failed  after  two  clearances, 
0.081  and  0.061  inches9  had  been  tested.  This  compressor  provides  air  to 
the  speed  controls  of  the  Allis-Chalmers  compressor  and  also  provides  the 
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air  which  is  used  as  a  reference  pressure  for  the  water  manometer  board. 
The  remaining  three  clearances  were  tested  using  a  smaller  auxiliary 
compressors,  which  introduced  considerable  fluctuations  into  the  house  air 
system. 

4.3  Theoretical  Analysis 

The  known  quantities  necessary  for  the  turbine  analysis  as  performed 
in  this  project  are;   mass  flow  rate;  total  inlet  temperature  and  pres- 
sure; static  inlet  pressure;  turbine  speed  and  torque  output;  atmospheric 
pressure  and  temperature;  and  the  absolute  rotor  inlet  flow  angle  and 
velocity  coefficient  based  on  the  results  of  the  scroll  and  guide  vane 
analysis. 

The  working  fluid  is  assumed  to  undergo  an  adiabatic  process,  since 
the  quantity  of  heat  passing  through  the  casing  is  negligible  compared  to 
the  enthalpy  change  undergone  by  the  fluid  as  it  passes  through  the  tur- 
bine.   It  is  additionally  assumed  that  the  fluid  behaves  as  a  perfect 
gas. 

The  flow  is  assumed  to  have  axisymmetric  stream  surfaces;  thus,  in 
effect,  the  analysis  is  one -dimensional,  based  on  the  determination  of 
the  flow  properties  along  a  particular  streamline  lying  on  a  chosen  sur- 
face of  revolution.   Vavra  (8)  states  that  the  major  portion  of  the  rotor 
losses  is  due  to  secondary  flow  phenomena,  tip  leakage  losses,  and  the 
loss  that  occurs  because  of  the  motion  of  the  blade  tips  in  the  wall 
boundary  layers.   Thus  is  would  appear  that  the  majority  of  the  losses 
occur  near  the  shrouds,  and  the  representative  streamline  was  accordingly 
chosen  as  that  existing  at  the  outer  radius  of  the  rotor  discharge.   It 
was  felt  that  this  approach  would  yield  a  better  correlation  of  results 
than  that  obtained  from  use  of  a  calculated  mass-averaged  mean  radius  at 

the  rotor  discharge  (7). 
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It  is  assumed  that  the  relative  flow  ani 


calculated  from  the 


lgle  ?    20  = 
blade  dimensions9  is  very  nearly  equal  to  the  actual  relative  flow  angle 

at  the  outer  discharge  radius  (where  R __  assumes  its  maximum  value).   This 

implies  that  the  so-called  slip  factor  may  be  neglected.   Verification  of 

the  validity  of  this  assumption  is  given  by  the  flow  surveys  taken  at  the 

rotor  discharge  by  Riley  (3).   Accordingly,  6  9n  was  calculated  from 

Equation  A2(4)  of  (1)  and  found  to  be  equal  to  -69.85  degrees  at  the  outer 

radius.   This  value  was  used  in  all  calculations. 

The  thermodynamic  process  of  a  fluid  passing  through  a  turbine  is 
given  in  Figure  14.   The  resulting  flow  velocities  are  shown  in  Figure  15. 
The  derivations  which  follow  are  based  on  these  two  diagrams. 

Referring  to  Figure  14,  for  a  reversible,  adiabatic  process, 


V    _   /   * 


y-i 


Tt 


?*, 


(16) 


where  T  ' '   represents  the  static  temperature  at  the  rotor  discharge  for 
an  isentropic  expansion  through  the  turbine  to  p..   The  isentropic  tempera- 
ture drop  which  would  occur  for  such  an  expansion  is 


aTls  =     Tro  -  Ta 


(17) 


and  may   therefore   be   expressed  as 


AJls    =  Tt< 


l-U/Pto)V 


(18) 


Similarly,    for  an   isentropic  expansion  across   the   inlet   guide  vanes 


Tt 


Pi. 


(19) 


and 


(AT«;) 


guide  vanes 


=  ho-Ji  =  Tii 


l-(Pl/p£o) 


Y-' 


(20) 
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The  isentropic,  or  theoretical  degree  of  reaction,  r,  represents 
that  portion  of  the  overall  isentropic  enthalpy  drop  which  occurs  in  the 
rotor,  or 

~T  '       T   " 

r    E       L  ~    z  (2D 

Therefore, 

r  -   )  —        s  &vi^e  vdn<5  (22) 

AT,S 
and,  from  Equations  (18)  and  (20), 

r  =  I  -    U~  (^/PiQ^J  (23) 

[l-  Vn/ru)*] 

The  total-to-static  efficiency  of  the  turbine  is  defined  as  the 
ratio  of  the  work  actually  produced  by  the  turbine  to  that  work  which 
could  ideally  be  accomplished  by  isentropic  expansion  from  total  inlet 
to  static  outlet  conditions „   For  this  case,  the  efficiency  may  be  ex- 
pressed by  the  ratio  of  the  corresponding  temperature  drops, 


Tl  =  ATw//aTC! 


(24) 


The  change  in  total  enthalpy  across  the  turbine  is 

m 

where  M  is  the  net  torque  produced,  1*J  is  the  rotational  speed,  and  W 
is  the  mass  flow  rate.   Therefore s    aTw  may  be  expressed  as: 

aTw  =   -^L_  (26) 

The  efficiency  given  by  Equation  (24)  was  computed  both  for  the 
case  where  bearing  friction  was  considered  and  for  the  case  where  it  was 
not.   Bearing  friction  was  calculated  based  on  the  results  of  coast-down 
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tests  performed  by  Vavra  (1).   The  friction  horsepower  is  expressed  by 

HPf  =  -  0.6  +  (N/iOOo)  (.17)43)  (27) 

for  N  <  10 9 500  rpm,  and  by 

HPf-  -0.6f  (hl/iOOo)Cnm)-(H.S98^ldiXH/lOOO-i0.5)2'    (28) 

for  N  >  10,500  rpm.   The  bearing  friction  moment  is  then  given  by 

Mr  =    HPf  '  55°  (29) 

The  net  moment  is  then  the  sum  of  T,  the  moment  measured  by  the  dynamo- 
meter, and  the  bearing  friction  moment: 

M  r  7  t  Mf  (30) 

The  velocity  C_  corresponding  to  the  isentropic  enthalpy  drop 
through  the  turbine  is  given  by  the  expression 

CoZ    5     ZgJcp   &TLs  (3i) 

and  the  peripheral  speed  of  the  rotor  at  the  rotor  inlet  is 

Uj  =   — — (32) 

1  720 

where  D..  is  the  rotor  inlet  diameter  in  inches .   By  Equations  (31)  and 
(32),  the  isentropic  head  coefficient  is  given  by 

The  absolute  velocity  V   is  expressed  by  Equation  (4).   An  alternate 
expression,  obtained  by  substitution  of  Equation  (23)  into  Equation  (4), 
is 
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Vi    =    <f>    7  ZgTcp   (l-r)  aTCs' 


(34) 


The  temperature  T  may  then  be  obtained  from 
and  M- ,  the  Mach  number  at  the  rotor  inlet,  is 


M,  =    )L 


i 


(36) 


Z$TcPTr 

The  results  of  the  dummy  rotor  tests  can  now  be  used  by  introducing 
Equation  (14)  to  obtain  ^  by  an  iterative  procedure  which  is  discussed 
in  Appendix  B.   Equation  (12)  gives  the  nozzle  loss  coefficient  as 

If n  =  1  -  4>  *■  <37> 

Values  for  cK±     can  now  be  taken  from  Equation  (13)  to  determine  the 
peripheral  and  meridional  velocity  components  of  V. . 

Vui   -  VI  Sin  oq  (38) 

Vmi   =  VL  COS  ©CA  (39) 

From  Figure  15  it  is  apparent  that  the  peripheral  component  of  the 
relative  rotor  inlet  velocity  VL  is 

Ww  =  Vui  -  u»  (40) 

and  W     is   then   found    from 

Wi      =  VVmi2     +   Ww*'  (41) 
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The  relative  inlet  flow  angle  K^  is 


pi     =     sin"1  (Wui/Vi) 


(42) 


At  this  point  all  conditions  at  the  entrance  to  the  rotor  are  known, 

2 
It  is  assumed  that  only  the  kinetic  energy  V   /2gJCp  is  useful  in 

the  rotor  (since  it  has  radial  blades),  and  that  the  incidence  loss  at 

2 
the  rotor  inlet  is  taken  as  W  -  /2gJCp  (8).   Thus,  the  effective  rela- 
tive velocity  head  is  reduced,  and  static  conditions  at  the  rotor  inlet 
are  really  those  at  T*  on  Figure  14. 

The  equivalent  state  point  E  shown  on  Figure  14  is  established  from 
the  energy  equation  for  a  rotating  system  in  the  following  manner  (6): 

H,l  =  ««W=  k  +  TgY  -  -jgy    =  hit  ^  -^f      <«> 

where  H  is  the  relative  total  enthalpy.  This  equation  may  be  rearranged 
R 

in  the  form 


He  rhi  +  ^pr  -  f   ^y  ^   -  h2 1  T^~ 


(44) 


where  the  left  side  is  equal  to  H  ,  the  total  enthalpy  at  the  state  point 

E 

E.   Equation  (44)  clearly  shows  that  state  point  E  is  the  proper  refer- 
ence point  from  which  to  view  the  expansion  process  through  the  rotor  in 
an  analogous  manner  to  that  of  the  stator.   The  peripheral  speed  at  the 
outer  discharge  radius,  U   »  is  given  by  Equation  (32)  with  D.  replaced 
by  D   ,  the  outer  diameter  of  the  rotor  discharge. 

The  work  output  of  the  turbine,  represented  by  the  change  in  total 
enthalpy,  may  be  represented  by  Euler's  turbine  equation  (6): 

4HW     =       C,ATW      r       Ul  VU1     -    U"  Vw  (45) 
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where  U„_  and  V  ~n  occur  at  the  outer  discharge  radius.   Manipulating 
this  equation  and  Equations  (24)  and  (31),  an  expression  for  the  peri- 
pheral component  of  the  absolute  discharge  velocity,  V TI9n»  is 


Vux< 


Rj 


vL 


R 


10 


Vtu n  /  _c 


u 


Ui 


(46) 


Using  the  velocity  diagram  of  Figure  15,  it  can  be  shown  that  the 
relative  velocity  at  the  outer  discharge  radius,  W_  ,  may  be  expressed 


by 


Wag  _.  (Vuzo  I  Vj)    -  (fa,/  Rx) 
Ui  Sin  (3: 


(47) 


2.0 


where,  in  keeping  with  the  stated  assumptions, 


£zo    =    -  69.  85 


(48) 


If  the  divergence  of  the  pressure  lines  on  the  entropy  diagram  of 
Figure  14,  in  the  direction  of  increasing  entropy,  is  neglected,  there  is 


T/-T.'  -   TV-  Ti 


(49) 


By  Equations  (21)  and  (31),  the  right  side  of  the  expression  is  equiva- 
lent to 


T  -  T    -    r 


(50) 


2-gTcp 
From  Figure  14,  the  left  side  may  be  expressed  in  terms  of  the  velocities 


T/-T* 


VJzot*    ,   UA*-1W 


V. 


mi 


(51) 


Z$Tcp  ^^JCp      Z%  Tcp 

Manipulation  of  Equations  (50),  (51),  and  (32)  yields  a  relation  for 


W„„  ,  ,  the  theoretical  relative  velocity  at  the  rotor  discharge: 
20th' 


w 


loth 


Vi 


IvjCoY 


m  ♦  m  - ' 


(52) 
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Using  Equations  (47)  and  (52)  it  is  now  possible  to  define  a  velocity 
coefficient  Y  f°r  tne  rotor,  which  is  analogous  to  the  velocity  co- 
efficient (f>    for  the  scroll  and  guide  vanes  given  by  Equation  (2). 

Y    =      Wzo/Wioih  (53) 

The  rotor  loss  coefficients  which  is  a  measure  of  the  kinetic  energy  loss 
through  the  rotor,  is,  then 


h  =  i  - .  y 


(54) 


The  meridional  velocity  component  of  the  absolute  discharge  velocity  can 
be  expressed  by 

Vmzo  -  Was  Zo5  fi,lo  (55) 

from  Figure  15.   The  absolute  discharge  flow  angle  is  then 

°^  2o   ~  tan  {yfuzo  /  \fmio  )  (56> 

If  consideration  is  given  to  the  case  where  bearing  losses  are  as- 
sumed to  exist,  the  resulting  efficiency  will  obviously  be  different 
from  that  for  no  bearing  losses.   Examination  of  the  foregoing  analysis 
shows  that,  correspondingly,  different  values  for  V  «,n>  ^M9n'  W?n8'  "^ 
and  Id      will  result. 

Thus,  from  the  measured  items  of  test  data  specified  at  the  begin- 
ning of  this  section,  the  following  performance  parameters  are  obtained 
by  the  analysis  technique;   flow  rate,  head  coefficient,  degree  of  re- 
action, absolute  and  relative  rotor  inlet  flow  angles,  velocity  and  loss 
coefficients  for  the  scroll  and  guide  vanes,  the  net  torque  produced,  the 
absolute  discharge  flow  angle,  the  rotor  velocity  and  loss  coefficients, 
and  the  total-to-static  efficiency. 
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In  this  experiment  the  above  parameters  are  referred  to  standard 
air  in  accordance  with  the  NASA  method,  as  described  in  (1).   This  method 
utilizes  a  total  inlet  pressure  of  14.7  pounds  per  square  inch  absolute, 
a  total  inlet  temperature  of  518.7  degrees  Rankine,  a  specific  heat  ratio 
of  1.4,  and  a  constant  pressure  specific  heat  of  0.24  BTU  per  pound  mass- 
degree  Fahrenheit.   The  referral  parameters  are 

I     =       Pro/  R.7  (57) 

9  E    Tto/sig.7  (58) 


and 


Yi 


^      _        Q.glO     /     Y4I\    2(V-Q  (59) 


The  factor  £  ,  which  corrects  the  flow  rate  through  the  turbine  for 
varying  values  of  y  ,  is  that  given  by  Vavra  (1)  and  is  supposed  to  give 
a  better  correlation  than  the  expression  for  €  used  by  NASA. 
The  referred  speed  is  then  given  by 


Nref.    =    N/fe  <60> 


The   referred   flow  rate   is 


W«f.   =  V/(/e"-f)  <61> 

Finally,     the   referred  moment   is 

Mref      E      M(-f)  (62) 
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4.4  Data  Reduction:  Programs  RADIAL  and  NOSPD 

The  computer  program  RADIAL,  listed  in  Appendix  B,  determines  the 
turbine  performance  parameters  described  in  Section  4.3  from  the  test  data. 
In  general,  the  specific  techniques  used  by  the  computer  in  these  computa- 
tions are  apparent  from  the  program  itself.   However,  Appendix  B  describes 
each  subroutine  briefly  and  elaborates  on  specific  computational  pro- 
cedures which  are  not  immediately  obvious. 

Program  NOSPD,  listed  in  Appendix  C,  was  written  to  compute  the  de- 
sired performance  parameters  from  the  locked  rotor  test  data. 

4.5  Discussion  of  Results 

The  input  data  and  results  of  program  RADIAL  are  listed  in  Tables  B4 
through  B13,  with  input  and  output  for  a  given  clearance  appearing  in  ad- 
jacent tables.   The  input  data  and  results  of  program  NOSPD  appear  in 
Table  C2. 

Figure  27  shows  the  degree  of  reaction,  r,  as  a  function  of  head  co- 
efficient, k.  ,  for  the  clearances  tested  in  this  analysis  as  well  as 
*   is'  J 

those  of  Riley  (3)  and  Vavra  (1).  For  each  clearance,  the  degree  of  re- 
action was  found  to  be  a  unique  function  of  head  coefficient,  independent 
of  turbine  pressure  ratio.   This  is  in  agreement  with  the  findings  of 
both  Riley  and  Vavra.   Riley  further  concluded  that  the  degree  of  re- 
action is  independent  of  axial  clearance.  The  results  of  the  present  tests 
verify  this  for  all  clearances  except  for  0.081  inch.  At  this  clearance, 
the  values  of  r  obtained  were  found  to  be  approximately  1  to  1.5  per  cent 
lower  than  those  for  the  other  clearances.   This  is  shown  on  Figure  27, 
where  the  resulting  curves  for  all  clearances  except  0.081  inch  are  rep- 
resented by  that  for  0.061  inch,  and  the  curve  for  0.081  inch  is  plotted 
separately.   This  discrepancy  implies  that,  for  a  given  overall  pressure 
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ratio,  the  pressure  p  can  no  longer  be  maintained  at  a  clearance  of 
0.081  inch. 

It  may  be  seen  from  the  representative  curve  of  Figure  27  that  r 
decreases  with  increasing  k,  (decreasing  speed),  from  a  value  of  ap- 
proximately  0.75  at  a  k.   of  1.05  to  0.30  at  a  k.   of  6.   This  drop  oc- 

IS  X  5 

curs  because,  with  decreasing  speed,  a  smaller  pressure  drop  is  required 

to  pass  the  flow  through  the  rotor. 

The  degree  of  reaction  curves  obtained  in  earlier  tests  by  Vavra  (1) 

and  Riley  (3)  are  superimposed  on  the  curve  of  Figure  27.   Vavra  states 

that  the  value  of  r  at  design  k„   of  2.3  should  be  approximately  0.486, 

as  compared  with  the  value  of  0.395  obtained  by  him  and  with  Riley's  value 

of  0.402.   The  current  tests  yielded  a  value  of  0.440  at  design  k.  , 

is 

except  for  a  clearance  of  0.081  inch,  for  which  the  corresponding  value 
of  r  was  0.434.   This  apparent  increase  in  degree  of  reaction  is  due  pri- 
marily to  a  more  accurate  determination  of  p.  than  was  used  by  either 
Vavra  or  Riley.   Boshoven  (4)  relocated  the  shroud  pressure  taps,  which 
measure  the  pressure  p.,  to  a  position  axially  outward  in  the  annulus 
around  the  periphery  of  the  shroud.   This  was  done  to  reduce  the  effects 
of  local  perturbations  on  the  pressure  readings  and  greatly  improved  the 
uniformity  of  the  pressure  distribution  around  the  shrouds.   This  im- 
provement in  the  accuracy  of  p1  yields  a  more  accurate  value  of  r  since 
degree  of  reaction  is  directly  dependent  on  p. . 

The  degree  of  reaction  remains  lower  than  the  design  value  because 
the  guide  vane  exit  area  is  too  small  to  permit  discharge  of  the  design 
mass  flow  (1). 

Significant  plots  resulting  from  the  overall  turbine  efficiency  cal- 
culations are  given  in  Figures  28  through  31.   Figure  28  depicts  turbine 
efficiency  as  a  function  of  referred  speed  for  the  case  of  no  bearing 
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losses  and  a  clearance  of  0.081  inch.   Figure  29  is  a  corresponding  plot 
of  efficiency  versus  head  coefficient.   Figures  30  and  31  represent  these 
same  respective  parameters^  also  at  a  clearance  of  0.081  inch,  for  the 
case  where  bearing  losses  are  considered. 

Efficiency  was  found  to  be  essentially  independent  of  clearance  over 
the  range  of  clearances  tested.  Figures  28  through  31  are  therefore  rep- 
resentative of  the  results  for  all  clearances  between  0.015  and  0.081  inch. 

Values  of  turbine  efficiency  were  obtained  from  measurements  of  the 
torque  at  the  turbine  shaft  made  with  the  air  dynamometer  and  from  the 
bearing  loss  values  as  determined  from  the  results  of  coast-down  tests 
conducted  by  Vavra  (1).   Because  of  the  previously  discussed  problem  of 
excessive  bearing  temperatures  and  seizing  experienced  during  these  tests, 
the  bearing  losses  are  probably  greater  than  the  maximum  losses  deter- 
mined by  Vavra.   However,  since  time  did  not  permit  a  repetition  of  the 
coast-down  tests,  Vavra' s  equations  for  maximum  bearing  losses  [Equa- 
tions (27)  and  (28)]  were  used. 

A  comparison  of  Figures  28  and  29  and  of  Figures  30  and  31  shows 
the  advantage  derived  from  representing  efficiency  as  a  function  of  head 
coefficient  instead  of  referred  speed.   Use  of  the  head  coefficient  ef- 
fectively reduces  a  family  of  curves  to  a  single  curve.   However,  if  the 
bearing  losses  are  ignored,  as  in  Figure  29,  the  resulting  plot  depends 
significantly  on  pressure  ratio.   In  this  figure  only  those  curves  rep- 
resenting the  pressure  ratios  for  maximum  and  minimum  efficiencies  are 
shown  to  avoid  clutter. 

Figure  31  shows  that  the  efficiencies  are  nearly  independent  of  pres- 
sure ratio  in  the  case  where  allowance  is  made  for  bearing  losses.   If 
these  bearing  losses  were  accurately  known,  it  would  appear  that  the  in- 
fluence of  pressure  ratio  would  be  completely  eliminated  from  a  plot  of 

efficiency  versus  head  coefficient. 
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In  general,  if  bearing  losses  are  not  considered,  a  pressure  ratio 

of  1.7  yields  a  maximum  efficiency  of  81.1  per  cent  at  a  k„   of  2.35  and 

is 

an  axial  clearance  of  0.015  inch.   Minimum  efficiencies  are  obtained  for 

a  pressure  ratio  of  1.2  and  are  in  most  cases  10  to  12  per  cent  lower 

than  maximum  efficiencies  in  the  vicinity  of  the  design  k,   of  2.3. 

is 

Vavra  (1)  reported  a  maximum  efficiency  of  81  per  cent  for  the  case  of  no 

bearing  losses  at  a  pressure  ratio  of  1.616.  k.   of  2.37 9  and  clearance 

is 

of  0.062  inch.   Riley  (3)  obtained  a  corresponding  maximum  efficiency  of 

82.52  per  cent  at  a  k„   of  2.26.  pressure  ratio  of  1.678,  and  clearance 

is 

of  0.027  inch. 

For  the  case  where  bearing  losses  are  taken  into  account,  the  maxi- 
mum efficiency,  ~f[L    ,  was  found  to  be  84.6  per  cent  at  a  pressure  ratio 
of  1.7,  head  coefficient  of  2.28,  and  clearance  of  0.081  inch.   This 
compares  favorably  with  the  design  efficiency  of  85.5  per  cent  at  a  pres- 
sure ratio  of  1.7  and  head  coefficient  of  2.3,  as  given  by  Vavra  (1). 
Riley  (3)  obtained  a  maximum  efficiency  of  86.04  per  cent  at  a  pressure 
ratio  of  1.678,  head  coefficient  of  2.259,  and  clearance  of  0.027  inch, 
based  on  the  maximum  values  of  bearing  losses. 

At  minimum  speed  (maximum  k„  )  for  a  given  pressure  ratio  and  clear- 
ance,  the  efficiency  in  general  is  about  65  per  cent.  As  speed  is  in- 
creased (k„   decreased),  the  efficiency  increases  to  a  maximum  of  about 

84  per  cent  in  the  vicinity  of  a  k„   of  2.3  and  then  falls  off  rapidly 

is 

to  30  to  40  per  cent  at  a  k„   of  about  1.05.   Initially,  then,  the  ef- 

X  s 

ficiency  increases  with  speed  because,  by  Equation  (26),  an  increase  in 
U  and  decrease  in  W  offset  the  decrease  in  torque,  resulting  in  in- 
creased turbine  work  and  efficiency.   As  speed  is  further  increased,  the 
decreasing  torque  begins  to  predominate,  causing  the  efficiency  curve  to 
peak  and  then  rapidly  decrease. 
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It  was  expected  that,  as  the  clearance  was  varied  over  the  range  en- 
compassed by  these  tests,  a  ""critical"  clearance  would  be  reached  beyond 
which  turbine  efficiency  would  deteriorate  rapidly.   Riley  (3)  observed 
this  to  occur  at  a  clearance  of  0.052  inch.   Such  an  occurrence  was  not 
observed  in  these  tests  over  the  range  of  clearances  from  0.015  to  0.081 
inch.   Thus,  it  appears  that  axial  clearance  has  no  appreciable  effect 
on  turbine  efficiency  over  this  range  of  clearances. 

It  has  been  previously  stated  that  Riley  (3)  was  in  all  likelihood 
operating  at  a  considerably  greater  clearance  than  he  believed  himself 
to  be.   This  is  substantiated  to  a  degree  by  the  fact  that  Vavra  (1)  did 
not  experience  this  efficiency  drop  at  a  clearance  of  0.062  inch.   Since 
gaps  also  existed  between  the  rotor  and  shrouds  during  Riley's  tests, 
the  blades  would  have  left  the  wall  boundary  layer  at  a  lesser  axial 
clearance  than  in  these  tests  where  the  blade  and  shroud  contours  were 
matched.   Thus,  it  is  reasonable  to  assume  that  the  critical  clearance 
is  at  some  value  greater  than  0.081  inch. 

That  such  a  drop  in  efficiency  should  occur  as  clearance  is  increased 
to  some  critical  value  is  explained  in  part  by  Csanady  (9).   The  rotating 
blades  "scrape  up"'  the  shroud  boundary  layer  and  produce  vortices  in  the 
vicinity  of  the  tips  of  the  blades  which,  in  turbines,  tend  to  nullify 
the  tip  vortices  themselves.   At  a  particular  tip  clearance,  then,  these 
two  vortices  would  neutralize  each  other,  leaving  only  the  passage  vortex 
which  causes  secondary  flow  losses.  This  neutralization  should  produce 
a  pronounced  increase  in  efficiency  at  that  clearance  at  which  it  occurs. 
For  the  subject  turbine,  it  would  appear  that  if  this  condition  does 
exist,  it  occurs  at  a  clearance  less  than  0.015  inch. 

The  efficiency  is  essentially  independent  of  clearance  because,  al- 
though the  tip  vortices  become  more  predominant  as  clearance  increases, 
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there  is  also  a  buildup  of  the  boundary  layer.   This  buildup  maintains 
a  reasonably  even  relative  balance  between  the  tip  vortices  and  the 
"scraped  up"  vortices.   As  clearance  is  extended  beyond  some  critical 
value,  however,  the  blades  begin  to  operate  near  the  edge  of  the  wall 
boundary  layer,  and  the  countering  vortices  quickly  cease  to  be  effec- 
tive in  neutralizing  the  increasing  tip  vortices.   The  result  is  a  pro- 
nounced increase  in  rotor  loss  and  a  corresponding  decrease  in  turbine 
efficiency.   It  is  also  worthy  of  note  that,  for  the  test  turbine,  al- 
though the  axial  tip  clearance  is  increased,  the  radial  clearance  between 
the  shrouds  and  the  rotor  blades  in  the  vicinity  of  the  discharge  remains 
essentially  unchanged. 

The  rotor  velocity  coefficient,  ^   ,  is  plotted  as  a  function  of 
the  head  coefficient  in  Figure  32,  for  the  case  where  bearing  losses  are 
neglected.   The  clearance  represented  is  0.024  inch.   Figures  33  through 

37  show  the  curves  of  Wl    versus  k   ,  with  bearing  losses  considered,  for 

is 

each  of  the  five  clearances  tested.  Again,  it  may  be  seen  that  if  bear- 
ing losses  are  ignored,  there  is  considerably  more  variation  of  the 
curves  with  pressure  ratio  than  in  the  case  where  bearing  losses  are 
considered. 

The  only  variation  in  rotor  velocity  coefficient  with  axial  clear- 
ance which  can  be  established  is  that  values  of  r   for  a  clearance  of 
0.081  inch  seem  to  be  generally  higher  than  for  all  other  clearances. 
This  variation  is  evident  for  both  bearing  loss  cases,  with  values  of 
^YL      being  generally  10  to  15  per  cent  higher  at  a  clearance  of  0.081 
inch.   Because  of  the  negligible  effect  of  axial  clearance  and  the  doubt- 
ful significance  of  the  case  where  bearing  losses  are  neglected,  Figure 
32,  for  a  clearance  of  0.024  inch,  is  the  only  curve  included  for  the 
case  where  bearing  losses  are  not  considered. 
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Figures  33  through  37  show  that  maximum  values  of  WL     do  not  always 
occur  at  the  higher  pressure  ratios,  although  Yl     is  more  nearly  con- 
stant over  the  range  of  k.   values  for  the  higher  pressure  ratios. 

X  s 

At  the  design  k.   of  2.3,  a  value  for  ^  of  0.960  was  produced  at 
is 

a  pressure  ratio  of  1.2  and  clearance  of  0.081  inch.  The  maximum  value 

of  "^   was  found  to  occur  at  the  same  pressure  ratio  and  clearance  and 

at  a  k„   of  3.35.   Its  value  is  0.990. 
is 

It  can  be  seen  from  the  curves  for  "%.      given  in  Figures  33  through 
37  that  Yl   remains  relatively  constant  over  a  range  of  k.   between 

about  2.3  and  the  maximum  value.  As  k.   is  decreased  below  2.3.   YL 

is  '   iL 

first  decreases,  with  a  minimum  occurring  at  a  k.   of  1.7  to  1.9.   It 

&       is 

then  increases  to  a  peak  at  a  k„   of  about  1.5,  and  falls  off  sharply 
as  the  head  coefficient  is  further  decreased.   Examination  of  the  tabu- 
lated performance  data  in  Appendix  B  shows  that  the  relative  rotor  inlet 

flow  angle  px  is  positive  at  values  of  k,   greater  than  about  2.3,  so 
1  is 

that  the  relative  flow  strikes  the  pressure  side  of  the  blades.  This 
produces  the  accompanying  incidence  losses,  and  further  losses  occur  as 
the  flow  is  turned  through  a  relatively  large  angle  to  follow  the  flow 
channel. 

As  IL  is  increased  (k„   decreased),  /Sj.  becomes  negative,  so  that  im 

1  is  ' 

pingement  of  the  flow  on  the  blades  acts  opposite  to  the  direction  of  ro- 
tation.  The  incidence  losses  are  not  significantly  changed;  however,  the 
flow  must  now  be  turned  by  a  considerably  smaller  amount  in  order  to 
follow  the  flow  channel.   Thus,  the  resultant  sum  of  these  three  losses 
appears  to  be  less,  and  an  increase  in  yL    occurs  at  a  k.   of  about  1.5. 

X  s 

Below  this  point,  as  Si  becomes  increasingly  negative,  the  relative  flow 
becomes  almost  perpendicular  to  the  blades.  Thus,  the  advantage  is  lost, 
and  the  rotor  efficiency  deteriorates  rapidly. 
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A  plot  of  referred  moment  with  bearing  losses  included,  as  a  function 
of  referred  speed 9  is  shown  in  Figure  38  for  data  obtained  at  a  clearance 
of  0.081  inch.   The  values  of  the  torque  at  zero  speed  resulting  from  the 
locked  rotor  tests  are  also  shown.   The  effects  of  clearance  on  this  para- 
meter were  also  found  to  be  negligible  over  the  range  of  clearances  from 
0.015  to  0.081  inch,  there  being  only  slight  random  variations  present  in 
the  resulting  data.   Therefore 9  Figure  38  is  a  representative  plot  for 
all  clearances. 

As  expected ,  the  maximum  moment  of  about  16.1  foot-pounds  was  found 
to  occur  at  zero  speed  and  at  a  pressure  ratio  of  1.7.   For  the  maximum 
speed  of  about  189000  rpm9  this  pressure  ratio  yielded  6.3  foot-pounds 
for  zero  bearing  losses  and  6„8  foot-pounds  for  maximum  bearing  losses. 
The  lowest  pressure  ratio9  1.2,  consistently  yielded  the  lowest  values 
of  moment.   The  moment  curves  of  Vavra  (1)  and  Riley  (3)  are  almost 
identical  to  those  obtained  in  these  tests,  except  for  the  results  obtained 
by  Riley  beyond  the  "'critical1"  clearance. 

In  general,  the  torque  produced  by  the  flow  energy  at  a  given  pressure 
ratio  is  a  maximum  at  zero  speed  and  decreases  with  increasing  turbine 
speed.   This  may  be  explained  by  the  law  of  angular  momentum,  which,  by 
reference  to  the  velocity  diagram  of  Figure  15,  can  be  written  as 


Al  =    VI  R»    [  (-1^)  Vm  -  V« 


(63) 


For  a  given  pressure  ratio,  V  .  is  nearly  constant  with  speed,  and  the 

* 

mass  flow  rate  W  decreases  with  increasing  speed.  From  Figure  15,  V  9 
also  increases  with  speed.  Thus,  it  is  seen  that  the  torque  produced, 
M,  decreases  with  increasing  turbine  speed. 
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After  completion  of  the  test  runs  previously  described9  an 
informal  trial  run  was  made  at  a  clearance  of  0.100  inch  in  an  attempt 
to  determine  the  critical  clearance.   Although  not  conclusive,  the 
results  indicated  no  significant  change  in  the  values  of  the  perform- 
ance parameters.   This  implies  that  a  reasonable  overall  efficiency  may 
be  maintained  at  clearances  which  might  considerably  exceed  those  tested 
in  this  analysis.   Since  these  additional  tests  were  not  complete,  the 
resulting  data  are  not  included  in  this  report. 
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SECTION  5 
CONCLUSIONS  AND  RECOMMENDATIONS 

At  a  pressure  ratio  of  1.7,  a  clearance  of  0.081  inch,  and  a  head 
coefficient  of  2.28,  a  maximum  total-to-static  overall  efficiency  of  84.6 
per  cent  was  obtained  for  the  case  where  bearing  losses  were  considered. 
The  efficiency,  referred  moment ,  and  degree  of  reaction  were  found  to  be 
independent  of  clearance  over  the  range  of  clearances  between  0.015  and 
0.081  inch.  A  clearance  tolerance  of  this  sort  is  very  important  in 
small  turbines  and  in  turbines  where  thermal  expansion  must  be  considered, 
It  is  recommended  that  additional  testing  be  performed  in  order  to  deter- 
mine the  precise  clearance  beyond  which  turbine  performance  deteriorates. 
Tests  might  also  be  conducted  in  which  the  clearance  at  the  rotor  dis- 
charge is  varied  in  addition  to  the  tip  clearance,  although  this  would  re- 
quire cutting  the  shrouds  or  the  manufacture  of  additional  shrouds  with 
various  annular  openings. 

The  results  of  this  investigation  show  that  computations  based  on 
flow  conditions  at  the  outer  discharge  radius  are  simpler  and  more 
straightforward  than  those  obtained  by  referring  the  losses  to  a  calcu- 
lated mass-average  radius,  and  that  this  simplicity  does  not  compromise 
the  quality  of  the  results.   On  the  contrary,  a  correlation  of  the  rotor 
velocity  coefficient  is  obtained  by  this  method,  whereas  previous  ef- 
forts to  obtain  meaningful  loss  parameters  were  not  successful  (7). 

There  is  strong  evidence  that  the  bearing  loss  data  used  in  this 
investigation  are  no  longer  applicable.   Correct  data  should  be  obtained 
either  through  new  coast-down  tests  or  through  some  implementation  of 
the  lubricating  oil  temperature  differential  across  the  bearings,  as  pre- 
viously described.   The  latter  method  would  be  highly  preferable  since  it 
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would  provide  accurate  loss  information  on  a  current  basis  and  would  cor- 
rectly reflect  any  change  in  the  mechanical  condition  of  the  bearings. 
Successful  implementation  of  such  a  method  would  probably  be  contingent 
on  a  complete  redesigning  of  the  test  rig,  so  that  the  oil  lines  are 
routed  outside  of  the  discharge  flow  as  much  as  possible  and  are  thoroughly 
insulated. 

Since  any  inaccuracies  in  torque  measurement  introduce  considerable 
error  into  the  calculations  as  performed  herein,  additional  effort  is  re- 
quired in  order  to  determine  how  to  properly  compensate  the  torque  capsule 
output  for  thermal  effects.  An  alternative  solution  would  be  to  utilize 
a  force  (reluctance)  capsule  and  a  moment  arm  arrangement.   This  capsule 
could  be  remotely  located  and  shielded  so  as  to  be  free  of  the  effects 
of  the  dynamometer  bearing  heat  and  the  cold  turbine  discharge  air. 

It  is  recommended  that  the  modifications  described  by  Vavra  (1)  be 
performed  in  order  to  improve  the  component  matching  of  the  turbine  and 
that  further  tests  be  carried  out  in  order  to  determine  the  effects  of 
these  modifications  on  the  performance  parameters.   Specifically,  the 
enlarging  of  the  guide  vane  discharge  area  or  the  reduction  of  the  rotor 
discharge  area  would  ultimately  increase  the  relative  rotor  discharge 
velocity  W-,  which  would  reduce  the  peripheral  velocity  V  _.   Greater 
turbine  work,  and  thus  higher  efficiencies,  should  result. 

It  is  felt  that  the  turbine  test  rig  is  now  in  much  better  operating 
condition  than  for  any  previous  tests.  The  instrumentation  has  been  proven; 
and,  for  the  first  time,  the  axial  clearances  may  be  accurately  set. 
Accordingly,  further  testing  to  verify  and  expand  on  the  performance  para- 
meters obtained  in  this  investigation  is  recommended. 
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FIGURE  5 
CROSS -SECTION  OF  TURBINE 


67 


69 


70 


72 


i    i)    ©  hi  ©  if 

I?    Jr-£-'^l|[J__         ^ 


Z 

o 

M 
H 

2 


z 


w 
z 
l-l 


|S 

M 

o 

M 

H 


u 


73 


FIGURE  12 
MERIAM  MERCURY  MICROMANOMETERS 
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FIGURE  13 
BROWN  POTENTIOMETER 
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FIGURE  14 
TEMPERATURE -ENTROPY  DIAGRAM  FOR  RADIAL  TURBINE 
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FIGURE  15 
VELOCITY  DIAGRAMS  FOR  RADIAL  TURBINE 
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FIGURE  16 
VELOCITY  COEFFICIENT  VS.  MACH  NUMBER 
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FIGURE  17 
ABSOLUTE  ROTOR  INLET  FLOW  ANGLE  VS.  MACH  NUMBER 
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FIGURE  18 
VELOCITY  COEFFICIENT  VS.  MACH  NUMBER 
FROM  DERIVED  EQUATION 
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FIGURE  19 
ABSOLUTE  ROTOR  INLET  FLOW  ANGLE  VS.  MACH  NUMBER 
FROM  DERIVED  EQUATION 
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FIGURE  22 
ELECTRONIC  SPEED  COUNTER 
AND  VIBRATION  ANALYZER 
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APPENDIX  A 

PROGRAM  SCROLL 

Program  SCROLL  computes  the  losses  in  the  scroll  and  guide  vanes  and 
the  absolute  rotor  inlet  flow  angle  from  torque,  mass  flow  rate,  pressure, 
and  temperature  data  obtained  in  the  dummy  rotor  tests   The  program  is 
essentially  that  used  by  Riley  (3)  and  by  Boshoven  (A),  with  several 
changes  as  specified  below.   It  can  process  any  number  of  runs,  with  a 
maximum  of  ten  sets  of  data  per  run   A  listing  of  program  SCROLL  is 
given  in  Table  Al.   The  input  and  output  values  which  constitute  the 
results  of  the  eight  sets  of  dummy  rotor  test  data  obtained  are  pre- 
sented in  Tables  A2  through  A9. 

Al   Main  Program 

The  value  of  the  specific  gravity  of  mercury,  G,  ,  at  room  tempera- 
ture t   ,  was  determined  from 
rm 

for  temperatures  between  32°F  and  113°F.   The  factor  for  converting 

2 
inches  Hg  to  lb/ft   is 

Cfc  -  ^*4388JL*  (  Gfcg  /[Z.5*\)  (A2) 

The  value  of  the  specific  gravity  of  water  at  room  temperature  is 

Gwo  ■  0.»«tt«  +•  (f.O*05"1»  mo*)tm  -  (iS$$ 186  JL(o*)   tn*1    (A3) 

The  specific  gravity  relations  were  obtained  from  the  tabulated  data  of 
Reference  10. .  These  equations  are  slightly  difftrent  from  those  used 
previously. 
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A2   Subroutine  TEMP 

The  relations  for  the  evaluation  of  the  temperatures  in  this  sub- 
routine from  the  thermocouple  readings  in  millivolts  are  different  from 
those  previously  used  because  of  the  substitution  of  iron-constantan 
thermocouples  for  the  chrome 1-alumel  type.   Thus,  for  iron-constantan 
thermocouples  producing  a  millivolt  output  V,  and  using  a  cold  junction 
reference  temperature,  t  „s  the  temperature  relations  are 

t  =   t.cj  +  0-m  +  3B.77  V  -    0.4518  \JZ  (A4) 

for      t  t  I00°F 

t  -    tl]  +  1/252  +  V\.U  V  -  0. 1555  V*  (A5) 

for       100'F  <  t  f   ZQO*F 

Equations    (A4)    and   (A5)  were   obtained   from  data   tabulated   in  Reference 

11. 

A3   Subroutine  FLOW 

The  evaluation  of  the  flow  rate  using  the  measured  data  from  the 
particular  flow  orifices  employed  in  the  radial  turbine  test  setup  was 
carried  out  in  accordance  with  a  method  developed  originally  by  Vavra  (1), 
after  the  techniques  given  by  Stearns,  et  al,   (5).   The  values  thus  deter- 
mined by  Vavra  applied  to  a  sharp-edged  orifice  of  2.800  inches  diameter. 
Before  this  series  of  tests  was  begun,  the  flow  orifice  diameter  was 
measured  and  found  to  be  2.7965  inches.  Additionally,  the  conversion 
factors  and  the  value  of  the  specific  gravity  of  mercury  used  herein  are 
slightly  different  than  those  previously  used,  the  latter  being  13.54  vice 
13.59.   The  equations  given  below  are  thus  those  of  Vavra  (1),  with  these 
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modifications  introduced.   Computations  were  made  for  the  vena  contracta 
taps  only,  although  flange  tap  data  were  available,  because  these  data 
give  a  more  accurate  flow  rate  (1). 

The  relation  for  the  flow  rate  is 


ft  =  C  <*  Yi  Fr   /  *^hivc  ' 


(A6) 


where : 


C        -  factor  dependent  on  orifice  diameter,  type  of  pressure 
taps,  and  dimensional  units 

<?C  -  area  multiplier  to  account  for  thermal  expansion  of 

orifice 

Y.        -  factor  which  accounts  for  compressibility  effects 

F         -  Reynolds  number  correction  factor 

r 

P.        -  absolute  pressure  at  upstream  tap 
Ah-       -  pressure  differential  across  orifice 

IVvj 

T,        -  temperature  ahead  of  orifice 
For  vena  contracta  taps,  C  ■  0.9057. 
For  a  stainless  steel  orifice: 

oC  =  1.0  +  (1.85*I0'S)  (T*- 532.)  (A?) 

Yi  =  1.0-    0.351  -SJSttL.  (A8) 

Also, 

Yr    =    1.0  f         °-  °°[215  (A9) 

where 

X    =    0.8151    -^  (Aio) 

Z   =    1-9  +   (2.*MlO-1)(T*-  b59.69)  <A11> 

where  Equation  (All)  is  valid  for  air  between  50°F  and  300°F. 
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The  resulting  flow  rate  equation  used  herein  is  then  given  by 

W  =  "354   c*  Yi  Yr     / 


?j.vc    a  hive 


Th 


(A12) 


where:  W  is   in   lb    /sec 


m 


&h  are   in  in.   Hg. 


T,    is   in  °R 
4 

and  where  c<  „  Y  and  F  are  as  given  previously. 


A4   Subroutine  PRESS 

It  has  been  previously  mentioned  that  a  regulated  reference  pressure 
was  used  on  the  96-inch  water  manometer  board.   This  pressure  must  there- 
fore be  accounted  for  in  computing  p„ ,  the  absolute  static  pressure  at 
the  rotor  inlet.   This  is  done  in  the  equation 


1  '  !!>.  SH 


where  the  quantity  (P    +H    ' 
M      J    x   atm   refl 


(A13) 

represents  the  total  reference  pressure 


applied  to  the  manometer  reservoir. 

A5   Method  of  Assembly  of  Input  Data 

The  items  of  input  data  necessary  in  performing  computations  with 
program  SCROLL  are  given  below  in  the  order  in  which  they  must  appear 
on  the  data  cards.   Specific  formats  for  the  numerical  entries  are  given 
in  the  program  listing. 


Card  No. 
1 
2 
3 


Number  of  runs  (clearances)  to  be  processed 

Number  of  points  (pressure  ratios)  in  run 

Measured  pressure  differential  across  flow 
orifice j,  from  vena  contracta  taps  (in.  Hg) 
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Card  No, 


FORTRAN 


Description 


PUVC 

P5P 

PATM 
HREF1 

HI 

SR 

TRM 

V4 

V5 

TCJ 

STARE 

Bl 


Measured  pressure  upstream  of  flow  measuring 
orifice 9  from  vena  contracta  taps  (in.  Hg 
gage) 

Measured  static  pressure  at  turbine  inlet 
(in.  Hg  gage) 

Barometric  pressure  (in.  Hg) 

Measured  reference  pressure  applied  to  HO 
manometer  board  (in.  Hg  gage) 

Average  static  pressure  ahead  of  rotor 
(in.  H20  gage) 

Torque  scale  reading  (lb.) 

Control  room  temperature  (°F) 

Thermocouple  reading  ahead  of  orifice  (mv.) 

Thermocouple  reading  at  turbine  inlet  (mv.) 

Cold  junction  temperature  (°F) 

Torque  scale  tare  (lb.) 

Perpendicular  distance  between  shroud  inner 
extremities  (in.) 


Cards  3  and  4  comprise  the  input  data  for  one  test  point  and  are 
repeated  for  a  total  of  NPTS  test  points  for  each  run.   Each  new  run 
begins  with  card  2. 
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APPENDIX  B 


PROGRAM  RADIAL 


Program  RADIAL  computes  the  referred  flow  rate,  head  coefficient, 
degree  of  reaction,  relative  and  absolute  rotor  inlet  flow  angles,  velocity 
and  loss  coefficients  for  the  scroll  and  guide  vanes,  referred  moment, 
absolute  discharge  flow  angle,  velocity  and  loss  coefficients  for  the  rotor, 
and  efficiencies  from  the  torque,  speed,  mass  flow  rate,  and  the  pressure 
and  temperature  data  obtained  during  the  active  rotor  tests.  The  calcula- 
tion technique  generally  follows  that  given  in  the  theory  of  Section  4.3. 
It  can  process  any  number  of  clearances  (sets)  and  any  number  of  pressure 
ratios  (runs)  per  clearance  with  a  maximum  of  ten  different  speed  settings 
(points)  per  run.   A  listing  of  program  RADIAL  is  given  in  Table  Bl. 

In  program  RADIAL,  use  was  made  of  temporary  memory  disc  storage  of 
the  input  data  in  order  to  avoid  the  awkward  procedure  of  subscripting  the 
input  variables.   This  procedure  permits  the  processing  and  printing  of 
the  input  data  on  a  point-by-point  basis  and  the  presentation  of  the  results 
in  a  compact,  tabular  form.  Appropriate  comments  appear  in  the  main  pro- 
gram at  the  points  where  this  procedure  was  used.   The  additional  control 
cards  required  to  implement  this  technique  and  to  obtain  necessary  addi- 
tional memory  cort  storage  space  are  shown  in  Table  B2.  As  the  table  indi- 
cates, these  cards  are  inserted  just  ahead  of  the  data  cards  in  the  as- 
sembled deck. 

The  torque  calibration  data  used  for  all  runs  are  given  in  Table  B3. 
The  input  and  output  values  resulting  from  the  five  sets  of  test  data  are 
presented  in  Tables  B4  through  B13. 
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Bl   Subroutines  TEMP  and  FLOW 

These  subroutines  perform  the  same  functions  as  in  program  SCROLL, 
namely  those  of  obtaining  temperatures  from  the  thermocouple  readings  and 
calculating  mass  flow  rate,  respectively.   Subroutine  TEMP  has  been  ex- 
panded to  include  calculation  of  the  temperatures  of  the  bearing  lubricat- 
ing oil  at  the  inlet  and  at  the  discharge  of  each  bearing.  However  for 
reasons  discussed  in  Section  4,  these  values  were  not  used  in  this  experi- 
ment. 

B2   Subroutine  PRESS 

Subroutine  PRESS  establishes  the  total-to-static  pressure  ratios 

across  the  scroll  and  guide  vanes  (P  /p,)  and  across  the  turbine  (P  / 

to  i  to 

p  )  and  the  ratio  p. /p  .  An  iteration  process  is  necessary  in  order  to 
determine  an  average  value  of  the  total  pressure  at  the  turbine  inlet, 
P  .   Three  relations  are  used  to  accomplish  this,  these  being  the  gas 


law, 


a      =  G,  — (bi) 


w< 


the  continuity  equation, 


V   -     — — (B2) 

and  the  energy  equation, 

where  A_  is  the  area  of  the  five-inch  inlet  pipe  and  C_  is  the  specific 
heat  of  air  at  the  temperature  T  .   Using  T  .  as  a  first  approximation 
of  P    ,  an  initial  value  of  V_  is  found,  and  from  this,  by  Equation  (B3) , 
a  new  value  of  T  .   The  iteration  continues  until  a  difference  of  0.01 


120 


degrees  or  less  exists  between  any  two  successive  values  of  T  .   Since 

velocities  are  low,  the  total  pressure  P  -  is  then  obtained  from 

tO 

B3   Subroutine  EDC 

This  subroutine  computes  the  mean  values  of  Y   and  C_  based  on  the 

arithmetic  mean  of  total  turbine  inlet  and  static  isentropic  discharge 

temperatures.   From  Reference  3,  the  variations  of  Y   and  C  with  tempera- 

P 

ture  are 


V=  1.4018  -  (2*  \0's)   t  <B5> 

CP  =  0.239H3  f  (3-4  no"')  t  4-  (2  *\CT§)   t1    (B6) 

These  relations  were  originally  obtained  from  data  tabulated  in  Refer- 
ence 12  and  are  valid  from  approximately  40°F  to  170°F. 

The  values  of  V    and  C,,.   v   correspond  to  the  average  tempera- 
0  avg      P(avg) 

ture  through  the  turbine  based  on  an  isentropic  expansion.   This  average 
temperature,  t,  is  given  by 


t  -  ("ft.  -  ^59.  <^) 


ATj, 


(B7) 


where 

r-i 


ATrt  =   Tti 


L-  (T*/Pt.) 


(B8) 


Using  an  iteration  process,,  the  first  approximation  of  4T,   is  based  on 
a  value  of  Y  corresponding  to  T  _.   The  value  of  t  resulting  from  Equa- 
tion (B6)  is  compared  with  the  previous  value  of  t  in  each  iteration 
until  the  difference  between  any  two  successive  values  of  t  is  less  than 
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1.0°.  Using  the  final  value  of  t,  values  of  y&^^  Cp(avg)s  and  dTis 
are  computed  from  Equations  (B5)  ,  (B6) ,  and  (B8)  ,  respectively.  The 
reference  parameter  9,  given  by 

©  =  (y«.,/M)  (Ts/ 5,8.7)  <B9> 

is  also  computed  in  this  subroutine. 

The  subscript  indicating  the  average  values  of  o  and  Cp  will  be 
dropped  in  subsequent  discussions.   The  average  values  of  these  quanti- 
ties are  assumed  in  the  remainder  of  the  program  unless  otherwise  stated. 

B4   Subroutines  TORQ  and  DYNA 

These  subroutines  compute  the  net  torque  produced  by  the  turbine  by 
evaluating  the  bearing  friction  moment  and  the  torque  measured  by  the 
dynamometer. 

In  subroutine  TORQ,  the  bearing  friction  moment  is  evaluated  from 
turbine  rpm  by  Equations  (27),  (28),  and  (29).   The  net  moment  is  calcu- 
lated by  Equation  (30),  where  the  dynamometer  torque  output  is  found  in 
subroutine  DYNA,  called  from  subroutine  TORQ. 

Subroutine  DYNA  determines  the  torque  from  the  dynamometer  calibra- 
tion data  obtained  for  each  run.   These  data  are  tabulated  in  Table  B3 
for  all  runs. 

The  values  of  the  torque  calibration  data  constitute  a  one-dimensional 
array  with  five  elements,  representing  the  torque  scale  reading  for  loads 
of  zero  to  400  inch-pounds,  respectively,  taken  at  100  inch-pound  intervals, 
Thus,  TCD(l)  represents  the  torque  scale  readout  value  for  the  no  load 
condition.   Using  a  DO  loop  with  index  I,  the  torque  indicator  reading 
for  the  test  point  being  processed,  TQ,  is  compared  with  the  values  of 
TCD(I)  to  determine  the  interval  (I,  1+1)   within  which  TQ  lies.   Since 
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the  calibration  curve  data  are  very  nearly  linear,  the  torque  is  computed 
assuming  a  straight  line  approximation  between  TCD(I)  and  TCD(I-l),  by 
the  equation 


T  -  ioo(i-i)  +  [oo 


TQ   -  TCDCr) 


TCBCr+1)-  TCD  (I) 
This  procedure  was  adapted  from  that  used  by  Riley  (3). 


(BIO) 


B5   Subroutine  TURB 

The  turbine  performance  parameters  described  in  Section  4  are  cal- 
culated by  subroutine  TURB  in  a  manner  which  follows  the  theory  previously 
given.   No  iterations  are  involved,  and  the  equations  used  are  self- 
explanatory.   Therefore,  the  steps  will  not  be  discussed  here.   The  veloc- 
ity coefficient,  (h    ,  and  the  absolute  rotor  inlet  flow  angle,  c<.,  are 
provided  from  subroutine  PHIAL,  for  which  TURB  is  the  calling  subroutine. 

B6   Subroutine  PHIAL 

This  subroutine  calculates  the  values  of  ^  ,  o(,9   V  ,  and  Lo     , 
from  the  results  of  program  SCROLL.   The  values  of  Y  and  Cp  used  are 
based  on  the  total  inlet  temperature,  rather  than  the  average  tempera- 
ture as  computed  in  subroutine  EDC,  because  the  values  of  p  and  o(  1 
obtained  in  program  SCROLL  were  computed  in  this  manner. 

The  values  of  (f>    and  oC     are  obtained  from  Equations  (13)  and  (14), 
for  which  the  entering  arguments  are  CLNC,  the  axial  clearance,  and  ACH1, 
or  M  ,  the  Mach  number  at  the  rotor  inlet.   Since  V-  and  T  are  not  known 
at  this  point  in  the  computations,  they  must  be  obtained  by  iteration. 
An  initial  value  of  V  is  obtained  from  Equation  (4)  by  assuming  (p  =  1.0, 
and  initial  values  of  T  and  M  are  obtained  from  Equations  (35)  and  (36), 
respectively.  Using  this  value  of  M  ,  called  ACHIA,  and  CLNC,  a  value  of 
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0,  called  PHIA,  is  obtained  from  Equation  (13).   The  assumed  value  of  0 
is  reduced  in  increments  of  0.001  until  succeeding  values  of  0  and  0 

A 

differ  by  no  more  than  0.001.   Thus,  values  of  V- ,  M  ,  and  0  are  estab- 
lished, and  o<  is  then  computed  from  Equation  (14).   The  loss  coeffi- 
cient for  the  scroll  and  guide  vanes,  w  ,  is  then  obtained  from  Equation 
(12). 

B7   Subroutine  REFER 

This  subroutine  computes  the  referred  speed,  referred  flow  rate, 
and  referred  moment  in  accordance  with  Equations  (60),  (61),  and  (62), 
respectively. 

B8   Method  of  Assembly  of  Input  Data 

In  order  to  facilitate  future  use  of  program  RADIAL,  the  order  of 
assembly  of  the  input  data  cards  is  given.   The  actual  formats  for  the 
numerical  entries  may  be  obtained  directly  from  the  program  listing. 
All  items  of  data  are  given  in  the  order  in  which  they  must  appear  on  the 
data  cards.   For  convenience,  the  symbols  are  redefined. 


Card  No, 


FORTRAN 


NSETS 


Description 


2 

CLNC 

3 

NRUNS 

4 

TCD 

5 

NPTS 

6 

RPM 

TQ 

PUVC 

Number  of  sets  of  data  (clearances  to  be 
processed 

Clearance  for  set 

Number  of  runs  (pressure  ratios)  in  data  set 

Torque  calibration  data  for  run 
(5  values) . 

Number  of  data  points  (speeds)  in  run 

Measured  turbine  speed  (rpm) 

Torque  indicator  reading  (counts) 

Measured  pressure  upstream  of  flow  measuring 
orifice,  from  vena  contracta  taps  (in.  Hg  gage) 
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Card  No.    FORTRAN 


P5P  Measured  static  pressure  at  turbine  inlet 

(in.  Hg  gage) 

HI  Average  static  pressure  ahead  of  rotor 

(in  H20  gage) 

DPVC         Measured  pressure  differential  across  flow 
orifice,  from  vena  contracta  taps  (in.  Hg) 

PATM         Barometric  pressure  (in.  Hg) 

HREF!        Measured  reference  pressure  applied  to  water 
manometer  board  (in.  Hg  gage) 

V4  Thermocouple  reading  ahead  of  orifice  (mv.) 

V5  Thermocouple  reading  at  turbine  inlet  (mv.) 

7       V20  Thermocouple  reading  at  lube  oil  inlet  (mv.) 

V21  Thermocouple  reading  at  lube  oil  discharge 

from  right  bearing  (mv,) 

V22  Thermocouple  reading  at  lub  oil  discharge 

from  left  bearing  (mv.) 

TRM  Control  room  temperature  (°F) 

TCJ  Cold  junction  temperature  (°F) 

Cards  6  and  7  comprise  the  input  data  for  one  test  point  and  are 

repeated  for  a  total  of  NPTS  test  points  for  each  run.   Each  new  run 

begins  with  card  4,  and  each  new  set  begins  with  card  2. 
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END 

/* 

//G0.FTO6FOO1    DD  SYSOUT=A,SPACE=(CYL, (15, 1)) 
//G0.FTO8FOO1   DD  UNIT=SYSDA,SPACE=(TRK, (2,1)) , 
//  DCB=(RECFM=V,LRECL=500,BLKSIZE=504) 

//G0.SYSIN  DD* 


16  COLUMN  72 


TABLE  B2 
ADDITIONAL  CONTROL  CARDS  REQUIRED  FOR  PROGRAM  RADIAL 
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TABLE  B3 
TORQUE  CALIBRATION  DATA 


SET 


CLEARANCE 
(IN) 


RUNS 


APPLIED  TORQUE  (IN-LB) 


0.081 


100 


200 


300 


1-4 

5-6, 
NOSPD  1 


5.    261, 


2.    255 


518. 


768 


400 


1013, 


510.    762.    1009 


0.061 


1-3 

4-6, 
NOSPD  2 


0.    255.     511.    760.    1007 


3.    259.     514.    762.    1007 


0.041 


1-3 

4-6, 
NOSPD  3 


0.    257 


511.    760.    1007 


0.    257.     510.    759.    1004, 


0.024 


1-3 

4-6, 
NOSPD  4 


0.    256.     510.    760.    1006, 


1.    259, 


513.    762.    1007 


0.015 


1-4 

5-6, 
NOSPD  5 


0.    255.     511.    760.    1008 


-2.    255. 


510.    759.    1004, 
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APPENDIX  C 


PROGRAM  NOSPD 


Program  NOSPD  calculates  the  referred  moments  and  referred  flow 
rates  for  the  stopped-rotor  tests  and  is  in  essence  a  much  simplified 
version  of  program  RADIAL.  The  computational  techniques  used  are  the 
same  in  both  programs  and  will  not  be  repeated.   Program  NOSPD  is  listed 
in  Table  CI,  and  the  input  and  output  date  are  given  in  Table  C2.   Table 
B3  lists  the  torque  calibration  data  used  for  each  run.   The  items  of 
input  data  defined  in  the  preceding  section  apply  to  this  program, 
except  as  noted  below. 

Description 

Number  of  runs  (clearances)  to  be  processed 


Card 

No. 

FORTRAN 

1 

NRUNS 

2 

CLNC 

3 

TCD 

4 

NPTS 

5 

V4,  V5, 
TCJ,  PATM, 
P5P,  DPVC, 
TQ,  TRM 
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